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1. INTRODUCTION 


BARYTES and celestite, sulphates of barium and strontium respectively, belong 
to the orthorhombic system with V,!° as the space group. The elastic con- 
stants of barytes have been determined by Voigt (1938) and by Bergmann 
(1938). Those of celestite have not so far been determined. In an attempt 
to study the elastic properties by the wedge method of single crystals belong- 
ing to the same isomorphous series, the elastic constants of the sulphates 
of barium and strontium have been determined by the author. 


2. THE ELASTIC CONSTANTS OF AN ORTHORHOMBIC CRYSTAL 


The orthorhombic system is characterised by nine constants C,,,Coo, 
Cy3, Cag, C55, Ces, Ciz, Cog and Cys. The a-, b-, c-crystallographic axes are 
respectively taken as the X,-, X.-, X,-axes. To get all the nine constants, 
at least six sections with different orientations are required. Six of the nine 
constants of the type C; come out directly from the axial cuts, using the 
secular Christoffel determinant. Thus a plate whose normal coincides with 
the a-axis (called the X,-cut) gives C,,, C;; and Cg; the X,-cut gives C2. 
C,, and C,, and the X;-cut yields Cz3, C,,and C;,; giving the six constants 
C;; with three checks. The remaining three constants Cy, C,,; and Cj; are 
determined by using suitable oblique cuts. For these cuts, however, a method 
which obviates the necessity of using approximate solutions is adopted. This 
method is just the same as the one given by Sundara Rao (1950) and hence 
only the essential equations are given below omitting the details. 


Taking the X, 45 X, cut (the cut whose normal is perpendicular to the 
X,-axis and is inclined at 45° to X, and X; and so on), the value of C4, is 
obtained from the relation 


C23 = + V {2x1 — (Cg3'+ Cys) P— 4 Cx”) —& 
Where x, is the longitudinal frequency of the particular section which 
can be experimentally determined. The rest of the quantities are previously 


determined, C3,’ +- Cy,’ being equal to 4 (Ca. + Cys + 2Cy,) and 4C;,’ 
=(C33— Cy»). 
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The ambiguity in the sign of the radical can, however, be decided easily 
since in the limiting case when C;,’-—»0 the value of C,3 should be that 
obtained by putting x, =Cg3,’, where C33’ = 4 (Cyo+ Cyg+ 2Co3+ 4C,4,) and 
whichever value of C,, is nearer the approximate one is taken as the correct 
one. 


In a similar manner, the other constants C,,; and C,, are determined 
from X,45 X; and X,;45 X, cuts using the relations given below: 


X2 45 X; cut: 
Cig= + V[f2x1— (Cos + Cy’)}* — 4Cg, 2] — C,; 
where Cy3'+ Cay’ = $ (Cy t+ Cagt+- 2C;5) 
and 4C 34’ = (Cgg3— Cy). 
In the limiting case, x, = C's; = 4 (C,,-+ Cgg+ 2Cy34- 4C,5). 
X; 45 X, cut: 
Cis= + VUf2%1— (Cas + Cyy’)}? — 4Cuu'*] — Coy 
where Cog’ + Cas’ = 4 (Cr+ Coot 2C ge) 
and 434° = (Cy — Cog). 
In the limiting case, x, = C33’ =4(Cy,+Co2+ 2Cyo+ 4C gg). 
3. Experimental Technique 


The crystal faces and axes are identified by using the crystallographic 
data given by Dana (1947). _ In the case of barytes the c-axis is easily identified 
being perpendicular to the prominent tabular face (001). The other two 
orthogonal axes a- and b- are located by examining the crystal between 
crossed polaroids. In the case of celestite, the faces are identified by measur- 
ing the interfacial angles with a contact goniometer. Both the crystals have 
an easy cleavage perpendicular to the c-axis and the X;-cuts are obtained 
by tapping the crystals with a sharp edge. 


All the six sections for each crystal are cut and ground to uniform thick- 
ness by hand and the different orientations are correct to within 1°. Wherever 
possible, more than one section of different thicknesses have been worked. 
The ultrasonic velocities in the crystal plates are determined by using the 
wedge method developed by Bhagavantam and Bhimasenachar (1944). 
Ultrasonic frequencies upto 12 Mc/Sec are used. The frequencies are deter- 
mined by a precision heterodyne wavemeter which is capable of giving an 
accuracy of | part in 50,000. 
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4. Experimental Data and Results 
The experimental data obtained and tne results calculated are given 
below: 
TABLE I 


Barium Sulphate i. Density = 4-432 gm./cm.® 
Unit : 10%! dynes/cm.? 

















i Pee al | Effective 
ro ye ona Thickness in mm. | frequency | Mode of elastic Value 
o. | of the plate | in Mc/Sec. vibration cman 
| 
i 4 X,-cot | 1-571 1-404 | i* Ci 8-622 
| do { do | 0-786 z - Cee 2- 703 
| do do | 0-812 t Css 2-885 
| | 
| 2 X,-cut 1-521 | 1-495 / Cus 9-166 
do do 0°822 t Cee 2-771 
do do | 0°550 t Cus 1-24) 
3 |. Xp-cut 1-509 | 1-639 ;| 10-84 
do do 0-842 t | Cs; 2-861 
| do | do | 0-535 t zo 1-155 
‘4 | X, 45 X2 | 1-689 | 1-259 1 | x, 8-017 
| & | X,45X5 | 1+872 1-236 i + 9-488 
' 
6 X, 45 Xy | 1-917 1-227 Z | mt 9-804 

















* Here and in the subsequent table, / and t denote the longitudinal and transverse modes 
respectively. 


TABLE II 


Strontium Sulphate Density = 3-955 gm./cm. 
Unit : 10% dynes/cm.? 











boi) | Fundamental | | Effective | 

|Plate | Description | > hickness in mm. Frequency Mode of | | elastic Value 

| No. | of the plate | ie Mc/See. | vibration oonmeent 

| 
eee Se < ee ee a Pn 
1 | Xi-cut | 1-766 | 1-455 ! Ci. | (10-44 
do do 0-758 t Css 2+834 
2 | X,-cut 2-061 | 257 l Sis 10-61 

do t Cee 2-726 


1 
do | 0-637 
do 0 
3 Xz-cut 1-766 
do 
do do 


4 | X,45X, 2-112 





2-194 


6 | X,45X; 





| 
| | 
| ° : vl | — ’ C66 2-585 
| | 
| 
2-625 | 1-040 | 








| 
| 
| 5 | X,45X; 
| 
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The values of the elastic constants calculated from the above data are 
given in Table ILI. 


Substance 


‘Barytes 


} 
| Celestite .. 





Ci 


8-62 
10-44 


TABLE III 
Unit: 10% dynes/cm.? 
| C22 Css Cas Cs; 
9°17 10-84 1-20 2-87 
10-61 12-86 1-35 2-79 











Cee 





2-74 


2+ 66 | 





Ci2 


5-23 
7-73 





Coz | Cys | 
3-56 | 3-41 | 
6-19 | 6-05 | 

| 





The elastic coefficients S’s are calculated using the conversion formule 
and are given in Table IV in units of 10-'* cm.?/dyne. 




















TABLE IV 

| | | | } 
| Substance Sit | S22 | S33 | Sas Ss5 Sec Siz S23 | S13 | 
| ® _ PEO, SOI SOR Tee Le = ss 
| | | | 
| Barytes ..| 18-40 | 17-36 | 10-96 | 83-33 | 34-84 | 36-50 —- 9-45|-2-73 |-2-68 

| 
|Celestite ..) 21-97 | 21-85 | 11-42 | 74-08 | 35-84 | 37-59 |—13-87|-4-00 |-3-66) 
| 





The linear compressibilities £,, 82, 8; along the principal directions 
X,, Xs, X; and the volume compressibility 8 in the above units are calcu- 
lated using the relations 


B;= Ba + Bj2 + Bs, where i=1, 2 or 3 and B=Z 8. 


TABLE V 
Unit: 10-15 cm.?/dyne 











| 
Substance Bx | Bo Bs B 
Barytes 6-27 5-18 5+55 17-00 
Celestite “| 4-44 3-98 3°76 12-18 





5. DISCUSSION OF RESULTS 


A comparison of the figures given in Tables III and IV shows that the 
values of elastic constants and elastic coefficients for celestite are higher than 
those for barytes. 


As has been mentioned already, the elastic constants of barytes have 
been determined by Voigt and Bergmann. Their values along with the 
author’s are given in Table VI. 
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TABLE VI 


Unit: 10%! dynes/cm.? 














Observer Ci Cos j Eas | C44 \ Css Coe Cis | Caz Cis 

Voigt --| 8°83 7°85 | 10-50 1-518 | 2-85 | 2-15 | 4-51 2-65 | 2-65 
| 

Bergmann ..] 8-83 | 7-81 | 10-38 | 1-17 | 2-80 | 2-55 | 4.77 | 2-80 | 2-69 
! 


Author ..| 8-62 | 9-17 | 10-84 | 1-20 | 2-87 | 2-74 | 5-23 | 3.56 | 3-4] 























An inspection of the above table shows that the author’s values are in 
general agreement with those of the previous workers but for C,. and the 
cross constants which are slightly higher. 


Table VII gives the values of the elastic coefficients S’s in units of 
10-13 cm.?/dyne calculated from the author’s C’s along with those of the 
other workers for comparison. The values of S’s of Voigt and Bergmann 
are the same as those quoted by Hearmon (1946). 


TABLE VII 


Unit: 10-3 cm.2/dyne. 





“55 66 


} ee ee ai a se 
Observer Si. | Soo | Ss Sa4 | S S S45 








: 
Voigt ..| 16-4 | 18-9 10-6 | 84-0 34-8 | 36-0 | -9-0 | -2-9 |-2-5 


Bergmann} 17-2 | 19-9 10-9 85+5 35°8 39-2 -9-9 | --1-7 |-3-0 











| 
| 
| | 
Author ..| 18°40 | 17-36 10-96 | 83-33 34 +84 36°50 | —9+45 | —2-68|—2-73 
| | 





The values of the compressibilities which have been computed from 
the above values of S’s along with the value of 8 obtained experimentally by 
Bridgman (1931) are given in Table VIII. 











TABLE VIII 
Observer Ba | Pe | Bs | p = =p | 
= : : + 
Voigt | 5:5 | ° Tod 6-2 Co 
Bergmann | 64 | 70 | G2 | 188 | 
Author | 6-27 | G18 | 5-55 | 17-00 | 
Bridgman oe oe ° | e | 17-63 | 
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It may be pointed out here that the author’s value is nearer the experi- 
mental value of Bridgman. 


The elastic constants of celestite have been determined for the first time. 
The value of compressibility for celestite determined by Bridgman (1928) at 
30° C. is available in the literature. The values of the principal compressi- 
bilities as computed from the author’s values of S’s lead to a volume 
compressibility of 12-18 x 10-'% cm.?/dyne as against Bridgman’s value of 
16-21 x 10-!%. This discrepancy cannot be accounted for. 


Table V shows that barium sulphate is more compressible than strontium 
sulphate. This may be due to the presence of the heavier barium ion in 
barytes. This is supported by the fact that the unit cell of barium sulphate is 
slightly larger than that of strontium sulphate. The dimensions of the edges 
of the unit cell are: 


Barium sulphate dy=8:89 A; by=S:'44A; Co=T:17A 
Strontium sulphate a)=8-36A; b)=5:36A; cy=6°84A 
The author takes this opportunity to offer his grateful thanks to 


Prof. S. Bhagavantam for his guidance and for the interest he has taken in the 
progress of the present work. 


6. SUMMARY 


The elastic constants of barium sulphate and strontium sulphate arfe 
determined using the wedge method. The values for barytes obtained by 
Bergmann and Voigt are verified. The elastic constants for strontium sulphate 
are determined for the first time and they are: 


C,, = 10-44; Cy. = 10-61; Cys = 12°86; Cy= 1-35; C,.= 2-79; 
Cog = 2°66; Cyg= 6°19; Cyg=6°05; and C,.= 7-73. 
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SYNTHETIC EXPERIMENTS IN THE 
BENZOPYRONE SERIES 


Part XVI. Isoformononetin and Its Derivatives 
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Received March 21, 1951 


IN two previous parts of this series! earlier experiments on the scope of 
partial methylation and demethylation for the preparation of 7-methyl ethers 
of isoflavones were reported with reference to the synthesis of prunetin and 
santal. The results were quite encouraging and both these naturally 
occurring 7-methyl ethers were prepared. In support of the theoretical 
considerations involved, it was felt that simpler cases should be studied in 
detail and the scope of the method fully understood. Hence the experiments 
described in this paper were undertaken. 


Of the two monomethyl ethers of daidzein (I), formononetin, the 
4’-methyl ether (II) was first found to occur in nature.2 This was first 
synthesised in poor yield by Wessely ef al.* by heating 2: 4-dihydroxy-4’- 
methoxy-phenyl-benzyl ketone with sodium and ethyl formate in a sealed 
tube at 100°. Subsequently Mahal, et a/.4 employed 2-hydroxy-4-benzyloxy- 
4'-methoxy-phenyl-benzyl ketone instead of the dihydroxy compound and 
carried: out the reaction at 0° and finally effected the removal of benzyl 
group. The protection of the 4-hydroxyl group in the ketone is not necessary 
and a good yield of formononetin has now been obtained directly from the 
dihydroxy ketone using the low temperature (0°) for the reaction. 
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The isomeric 7-monomethyl ether of daidzein which is now named 
isoformononetin (III) has not so far been made or described. For its 
preparation we have employed the partial methylation of daidzein using 
one mole of dimethyl sulphate. The product is a mixture which could be 
separated into the following fractions: (1) carbonate soluble; (2) carbonate 
insoluble but sodium hydroxide soluble; (3) sodium hydroxide insoluble. 
Fraction (1) obtained by acidifying the carbonate solution consisted of pure 
daidzein and no formononetin seemed to be present in it. Fraction (2) 
precipitated from sodium hydroxide solution by acid yielded isoformononetin 
and fraction (3) which was rather small yielded dimethyl daidzein. The 


results are in accordance with the theoretical considerations discussed in 
Part XMPII.! 


The method of partial demethylation was next investigated using O- 
dimethyl daidzein and hydrobromic acid and hydriodic acid; under suitable 
conditions almost quantitative yields of isoformononetin could be obtained 
with these reagents. Heating with hydriodic acid at 120° for half an hour 
seems to offer the most satisfactory method of partial demethylation for 
obtaining the 7-methyl ethers of isoflavones. 


5: 7-Dihydroxy and dimethoxy isoflavones offer the simplest examples 
for the study of partial methylation and demethylation. Based on analogy 
with similar flavones and flavanones the existence of large difference between 
the two hydroxyl groups could be readily understood as also the possibility 
of partial methylation of the 7-hydroxyl group. For obtaining the 7-methyl 
ether by demethylation, mild demethylating agents would normally be 
suggested just as in the case of chrysin dimethyl ether. But in view of our 
experience with other isoflavones we have now carried out this partial 
demethylation using hydriodic acid and the method works quite satis- 
factorily. 


5: 7-Dihydroxy isoflavone and its methyl ethers were not described 
till the present work had been completed and the recent paper of Iyer et al.® 
containing this description appeared only in March this year. Since our 
methods and results are slightly different, they are reported in this paper. 


2-Hydroxy-4: 6-dimethoxy-phenyl-benzyl ketone is conveniently made 
from phloroglucinol trimethyl ether and phenyl acetyl chloride by condensa- 
tion in presence of anhydrous aluminium chloride. The product agrees in 
all its properties with that described by Badcock, et al.* and Iyer, et al.’ who 
have made it by a different route. The subsequent isoflavone condensation 
proceeds quite smoothly giving rise to high yields of 5: 7-dimethoxy iso- 
flavone, Partial demethylation of this dimethyl ether with hydriodic acid 
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gives a very high yield of the 7-monomethyl ether. The interesting point 
that should be mentioned here is that in the analogous case of chrysin- 
dimethyl ether using similar conditions the product is found to be mostly 
chrysin. This reveals the existence of marked difference in the stability of 
the 7-methyl ether groups in flavones and isoflavones. The explanation is 
to be found in the difference in the location of the side phenyl nucleus in the 
two groups of compounds. In flavones the demands of the pyrone C-=O 
group for electrons are met by both the phenyl nuclei whereas in isoflavones 
the side phenyl plays no part. Consequently in the latter, the 7-position 
receives a higher quantum of positive charge making this particular hydroxyl 
more strongly acidic and the methoxyl more firmly held. 


EXPERIMENTAL 


Formononetin.—Powdered sodium (1 g.) was cooled in ice and treated 
with a solution of 2: 4-dihydroxy-4’-methoxy-phenyl-benzyl ketone (2 g.) 
in ethyl formate (10c.c.) with shaking. After keeping in the refrigerator 
for 48 hours pieces of ice were added and the excess of ethyl formate was 
removed under reduced pressure. On acidifying the cooled solution 
formononetin separated as colourless solid and crystallised from alcohol 
as colourless prisms melting at 255-57°.? It dissolved readily in aqueous 


sodium carbonate solution and gave no colour with ferric chloride in 
alcoholic solution. 


Daidzein—This was obtained by demethylating formononetin with 
hydriodic acid (d. 1-7) at 140° for one hour. Jt crystallised from alcohol 
as colourless prisms melting at 320-22°. The diacetate crystallised from 
ethyl acetate as clusters of needles melting at 184-6°.4 


Partial methylation (7-Methoxy-4'-hydroxy isoflavone, isoformononetin).— 
To a solution of daidzein (0-9 g.) in acetone (100c.c.) dimethyl sulphate 
(0:35c.c.) and potassium carbonate (1 g.) were added and the mixture 
refluxed for three hours. Acetone was then distilled off and the residue 
treated with water. The solid that separated was filtered and macerated 
with sodium carbonate. The residue was then treated with sodium 
hydroxide solution and the insoluble dimethyl ether was separated by 
filtration. The dimethyl ether (50 mg.) crystallised from alcohol as colourless 
needles melting at 154-56°. The alkali solutron on acidification gave a 
colourless solid which crystallised from alcohol as long stout elongated 
rectangular prisms melting at 218-20°. It was insoluble in cold sodium 
carbonate solution and gave no colour with ferric chloride in alcoholic 
solution. Yield, 0-4g. (Found: C, 71:2; H, 4:5; (Cj sH,,O, requires 
C, 71-6; H, 4°5%.) 
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The carbonate solution on acidification gave a colourless solid melting 


above 310° and one crystallisation from alcohol gave pure daidzein, m.p. 
320.-22°. 


Acetate.—7-Methoxy-4’-hydroxy isoflavone obtained above was acetylated 
by refluxing with acetic anhydride and pyridine for one hour. The 
acetate crystallised from ethyl acetate as colourless plates and prismatic 
needles melting at 148-49° (Found: C, 69-4; H, 4:4; C,,H,,0; requires 
C, 69-7; and H, 4-5°,.) 


On deacetylation with alcoholic hydrochloric acid and crystallisation 
from alcohol the isoflavone melted at 218-20°. 


7: 4’-Dimethoxy isoflavone.—2 : 4-Dihydroxy-4’-methoxy-phenyl-benzy] 
ketone required for this purpose was prepared following the method of Baker 
and Eastwood,’ using resorcinol and p-methoxy-benzyl cyanide. This was 
methylated to 2-hydroxy-4: 4’-dimethoxy-phenyl-benzyl ketone using dimethyl 
sulphate (1 mole) and anhydrous potassium carbonate in acetone solution. 
It (1 g.) was dissolved in ethyl formate (6c.c.) and the solution added to 
pulverised sodium (0-5 g.), cooled in ice and the mixture left in the refrigerator 
for 48 hours. Crushed ice was then added and ethyl formate removed under 
reduced pressure. On leaving overnight a colourless solid separated out. 
It was filtered and the filtrate acidified when some more of the product was 
obtained. 7: 4’-Dimethoxy isoflavone crystallised from alcohol as colourless 
needles melting at 154-56°, was insoluble in aqueous sodium hydroxide and 
gave no colour with ferric chloride. 


Partial demethylation 


(i) With hydrobromic acid.—7:4'-Dimethoxy isoflavone (0-3 g.) was 
dissolved in glacial acetic acid (5c.c.) and heated with hydrobromic acid 
(d. 1-51; 10 c.c.) for 2 hours on a boiling water-bath. The mixture was 
then diluted with ice-water and the solid that separated was filtered and 
washed with water. It was then dissolved in 5% aqueous sodium hydroxide, 
the solution filtered and the cold alkaline solution acidified with hydro- 
chloric acid. On crystallising from alcohol the product was obtained as 
colourless long rectangular prisms melting at 218-20°. It was identical 
with 7-methoxy-4’-hydroxy isoflavone in its properties and a mixed melting 
point with the sample described above was undepressed. 


(ii) With hydriodic acid.—7: 4’-Dimethoxy isoflavone (0:4 g.) was dis- 
solved in acetic anhydride (4c.c.) and treated with cooling with hydriodic 
acid (10c.c.; d. 1-7). The solution was heated in an oil-bath at 120° for 
half an hour, It was then diluted with water, iodine decolourised with 
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sulphurous acid and left in the refrigerator overnight. The solid that 
separated was filtered, washed with water and macerated with sodium carbo- 
nate solution. The residue was then dissolved in aqueous sodium hydroxide 
and the clear solution acidified. The product crystallised from alcohol as 
colourless stout rectangular prisms melting at 218° and was identical with 
7-methoxy-4’-hydroxy isoflavone; the mixed melting point with the above 
two samples was undepressed. Yield, 0-25 g. The acetate crystallised from 
ethyl acetate as colourless needles and melted at 148-49° alone or in admixture 
with the other samples already described. The sodium carbonate extract on 
acidification gave a very small quantity of daidzein. 


7-Methoxy-4'-ethoxy isoflavone.—7-Methoxy-4’-hydroxy isoflavone 
(0-1 g.) was refluxed in acetone solution with excess of ethyl iodide and 
anhydrous potassium carbonate for 6 hours. The ethyl ether crystallised 
from dilute alcohol as colourless short stout prisms melting at 121-22°. It 
was insoluble in aqueous alkali and gave no colour with ferric chloride in 
alcoholic solution (Found: C, 73:4; H. 5:2; CysH,;,O, requires C, 73-0; 
H, 5-4%.) 


Ethylation of formononetin (7-ethoxy-4'-methoxy isoflavone) 


Formononetin (0-1 g.) was ethylated as described above. Ethyl for- 
mononetin crystallised from alcohol as colourless stout rectangular prisms 
melting at 141-42° (Found: C, 72:8; H, 5:6; C,sH,,O,4 requires C, 73-0; 
H, 5-4%.) 

2-Hydroxy-4 : 6-dimethoxy-phenyl-benzyl ketone.—Phloroglucinol trimethy] 
ether required for this purpose was conveniently obtained by the dry 
methylation of phloroglucinol with dimethyl sulphate (3 moles) and anhydrous 
potassium carbonate in acetone solution. 


To a solution of anhydrous aluminium chloride (12 g.) and phloro- 
glucinol trimethyl ether (5 g.) in dry ether (100c.c.) cooled in ice, phenyl 
acetyl chloride (4.c.) was added drop by drop with shaking. The mixture 
was allowed to stand in the ice-bath for half an hour and then at room 
temperature for 24 hours. Ether was distilled off and the residue treated 
with ice and concentrated hydrochloric acid. After warming on a water- 
bath for 15 minutes the mixture was cooled. The oily liquid soon solidified. 
It was filtered and washed with dilute hydrochloric acid. The residue was 
dissolved in aqueous sodium hydroxide and when the solution was acidified 
it gave a colourless solid. On crystallisation from alcohol it separated as 
colourless prisms and melted at 117-18°.5° It gave a reddish brown colour 
with ferric chloride and a blue colour with concentrated nitric acid. 
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5: 7-Dimethoxy isoflavone——The above ketone (1 g.) was dissolved in 
ethyl formate (5c.c.) and cooled in ice. Powdered sodium (0-5 g.) cooled 
in ice, was then added in small quantities and the mixture left in the refrige- 
rator for 48 hours. Crushed ice was added and excess ethyl formate removed 
under reduced pressure. On cooling, the oily liquid soon solidified. It 
crystallised from alcohol as small prisms and melted at 121-2°. Some more 
of it could be obtained by acidifying the aqueous solution. Yield 0-6 g. 
Iyer, et al. gave its melting point as 120°. 


5: 7-Dihydroxy isoflavone—This was obtained by demethylating the 
dimethoxy compound with hydriodic acid for 3 hours at 140° in an oil-bath. 
5 : 7-Dihydroxy isoflavone crystallised from alcohol as colourless prisms 
melting at 192-4°. It readily dissolved in aqueous sodium carbonate and 


gave a pink colour with ferric chloride. Its properties agreed with the descrip- 
tion of Iyer, et al.5 


5-Hydroxy-T-methoxy isoflavone 


(1) Partial methylation—On methylating with one mole of dimethyl 
sulphate and anhydrous potassium carbonate in acetone solution for 3 hours 
5: 7-dihydroxy isoflavone gave the 7-methyl ether. It crystallised from ethyl 
acetate as colourless rectangular prisms melting at 140-41°. It gave a red 


colour with ferric chloride and was sparingly soluble in aqueous sodium 
hydroxide. 


(2) Partial demethylation.—S : 7-Dimethoxy isoflavone (0-4 g.) was heated 
with acetic anhydride (5 c.c.) and hydriodic acid (10c.c.) at 120° for half an 
hour and the product worked up as usual. It was macerated with sodium © 
carbonate solution and the residue crystallised. It separated from ethyl 
acetate in the form of colourless stout rectangular prisms and melted at 141°. 
(Found: C, 71:2; H, 4:3; CysH,2O, requires C, 71-6; H, 4:5%) It 
agreed in all its properties with the sample obtained by method (1) and a mixed 
melting point was undepressed. The carbonate washings when acidified did 
not yield any 5: 7-dihydroxy isoflavone. 


The above samples agreed with the description of Iyer et al.5 


SUMMARY 


The study of daidzein derivatives provides further support to the con- 
clusions already arrived at regarding the special features of partial methyla- 
tion and demethylation in the isoflavone series. Partial methylation of 
daidzein gives a good yield of the 7-monomethyl ether (isoformononetin) 
and the formation of the isomeric formononetin could not be detected. The 
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Same 7-methyl ether is obtained most conveniently by the demethylation of 
daidzein dimethyl ether with hydrobromic acid and hydriodic acid. 5:7- 
Dihydroxy and di-methoxy isoflavones have also been examined with a view 
to compare them with corresponding flavones. 7-Methoxy-5-hydroxy iso- 
flavone is obtained most readily by partial demethylation of dimethoxy iso- 
flavone with hydriodic acid whereas chrysin dimethyl ether yields mostly 
chrysin under these conditions. 
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1. INTRODUCTION 


No non-static solutions, with physical significance, of Einstein’s field equa- 
tions are known outside the field of Cosmology. The field of a radiating 
mass presents a problem for which general relativity has, so far, not been 
able to provide a solution. Schwarzschild’s external solution deals with the 
gravitational field of a cold dark body whose mass is constant. The applica- 
tion of this solution to describe the sun’s gravitational field should only be 
regarded as approximate. Various attempts have been made to generalize 
Schwarzschild’s solution in order to make it applicable to non-static masses, 
(Narlikar, 1936; Narlikar and Moghe, 1936). 


While discussing this outstanding unsolved problem of general rela- 
tivity, Professor Narlikar (1939) remarks: 


‘** If the principle of energy is to hold good, that is, if the combined 
energy of the matter and field is to be conserved, the system must be 
an isolated system surrounded by flat space-time. A spherical radiat- 
ing mass would probably be surrounded by a finite and non-static 
envelope of radiation with radial symmetry. This would be surrounded 
by a radial field of gravitational energy becoming weaker and weaker 
as it runs away from the central body until at last the field is flat at 
infinity. It has yct to be seen whether and how this view of the distri- 
bution of energy is substantiated by the field equations of relativity.” 


We represent below the solution of the field equations which substantiates 
the views expressed above. We begin with the derivation of the energy 
tensor for the radiation envelope surrounding a star. 


* The treatment as given here is essentially different from that of Professor H. Mineur as 
it appears in Ann. de l’Ecole Normale Superieure, Ser. 3, 5, 1, 1933. Our attention was 
kindly drawn to it by Professor Mineur some years ago.— V.V. N. 23-4-1950. 
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2. ENERGY TENSOR FOR A DIRECTED FLOW OF RADIATION 


By the term “directed flow of radiation’? we mean a distribution or 
electro-magnetic energy such that a local observer at any point of the region 
of space under consideration finds one and only one direction in which the 
radiant energy is flowing at the point. Using natural co-ordinates at the 
point of interest, we may take the components of the energy tensor as being 
given in terms of electric and magnetic ficld strengths E and H by the typical 
examples given by Tolman (1934). 


T,!? = 2 4 (E,?— E,*— E?+ H,2— H,?— H.?), (2-1) 
T,!? = — (E,E, + H,H,), (2-2) 
Ty4= (E,H.— E,H,), (2-3) 
Te“ = 4 (62+ E,'+ E,?+ H,*+- H,?+ H,’). (2-4) 


The suffix 0 to a component of a tensor indicates that the component is 
evaluated in natural co-ordinates at the point of interest. Considering, for 
simplicity, that the axes of our natural co-ordinates are oriented in such a 
way that the flow of radiation at the point of interest is in the x-direction 
and further that the radiation is polarised with the electric vector parallel 
to y-direction, we shall find 


E,=E,=H,=H,=0; E,=H., (2-5) 
and so the only surviving components of the tensor T,“” would be 
TU = T= T,4=4 (E,?+ H.”) =p, (2:6) 


p being the density of the radiant energy at the point. 


Having obtained the components of T+” for one system of co-ordinates, 
we can find them in any other system by the rules of tensor transformation. 
For a general co-ordinate system with the line-element 


ds* = g,,dx"dx", (2-7) 


the components of T#” will be given by 


wa w& . 
= St, xk, TF. (2-8) 


On using (2-6) this yields 
TH — [= ax”, (Ua G2” ax# Ox” =, ax dx” | P 


wv 


ke a oe og a ae 
IXy! Xp! «OX * OX" IX? IXp* | - AX op" AXgt 


(2-9) 


As the radiant energy travels along null-geodesics 
dx,' = dxo* = dr (say). 
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By (2-10) along the radiation flow we find 


Syydx*dx” = 0, (2-11) 
Next we use (2-10) in 

dx" dx dx,% 

dr _ dX" dr 
and find 

dx*  dx#  dx# 

dr dx! * dxoH° (2-12) 
With the help of (2-12), (2-9) finally reduces to 

1x” dx” 

Te = pm = : (2-13) 

with 
dx* dx” 
Suv dz dt aes (2:14) 


Thus for our case of the outside field of a non-static mass the energy 
tensor is to be taken of the form 


TH = putty”, (2-15) 
with 
14 =0; (0), 0” =0. (2-16) 


3. THE FIELD EQUATIONS 


A star of mass M and radius ry is supposed to start radiating at time fp. 
As the star continues to radiate the zone of radiation increases in thickness, 
its outer surface at a later instant ¢, beingr=r, Forrg<r<n,ty<t<t, 
let the line-clement be assumed to be of the form 


ds* = — e\dr?— r* (d6? + sin? Odd?) +- e’dt?, 


A=A(r, 0, v=v(r, t). (3-1) 
For the nature of radiation we have found the energy tensor T“” of the form 
"TH = pute’, (3-2) 
p is the density of radiation and the lines of flow are null-geodesics: 
v0 =0; (v4), 0” =0. (3-3) 
Since (T#”), = 0, we have the analogue of the equation of continuity 
(pv“),, = 0. (3-4) 


As the flow is to be radial, v?=0, v?=0 and 


7; = pv,v', T,f = pu sv, Ty — pv,v4, T,* => T,? = 0. (3 ° 5) 
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Also v,7“ =0 simplifies to 


—_ er (v})2 4. @” (v4)? =(). (3-6) 


With the usual expression for the components of T“” in terms of g,, and 
their derivatives, (3-5) gives the following three field equations: 


(i) T,tetr-> 24.7 4=0, 3 . 7) 
| v1 .. £ 

rv — - a Se o—(A4+-v)/ 2 — . J 
. ME ape pt re 6-8) 
(ii) T,'+- T,4=0, (3-9) 

V—v 2 2 

eS 5 See Se 4 4. . 
or (7 *) + 5=0: (3-10) 
(iii) T.2=0, (3-11) 

ap ¢ Cn Fo ht BM 
or e (5 ie a -) +¢ G ‘7 ~ eee (3-12) 


Here and in what follows an overhead dash or dot indicates a differentiation 
with regard to r or f. 


If the total energy is to be conserved, the line-element obtained by 


solving these equations must reduce to the static form 


=a 
pe (1 so SY dr? — +? (d02 +- sin? 0d4®) -- (1 i ) dt? (3-13) 


> 
at r=r,, t= and forr>r, at t =t, 
4. THE SOLUTION OF THE FIELD EQUATIONS 


2m 


On putting ec =1- , > m=m(r, bt) (4:1) 


in the field equation (3-8) we find that it is equivalent to 


ena SM 4 rit M0, (4-2) 
Using the operator 

iy t td (4-3) 
we may express this as 

dm _ 0. (4-4) 


dr 
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From (4:2) we can express e”? in terms of m: 


etn — = ( = omy" 2 (4-5) 


I 


Now we can take the second field equation (3-10). On substituting the values 
of A and v from (4-1) and (4-5), we find that 


nr m 2m on 
G 2 4 =5 pet (4-6) 
The first integral of the above equation is 
il (1 = =) — f(m), (4-7) 


f(m) being an arbitrary function. (4-7) is the differential equation to be 
solved for m. 


We now take the third field equation (3-12). We shall show that when 
A and v are given by (4-1), (4:5) together with the last differential equation 
(4-7), the equation (3-12) is automatically satisfied. The following is an 
identity holding between the components of the tensor Tes 


v’ 2 , Aa 
+ 9 - 5 eT) + F-  + TA(AF)= 0. 
(4-8) 


With the help of this identity and the two equations (3-7) and (3-9) the 
equation (3-11) can be transformed into 


J (PeA Ti) =0, (4-9) 


Thus the third field equation is satisfied, i.c., T,2=0 provided (4-9) is satis- 
fied, i.e., provided 

dj. 2m\) __ 

Ir yn (1 — )j =(, (4-10) 


i.e., provided “g = 0 when we use (4:7). And the last relation is already 


proved as (4°4) above. 


Hence we have solved all the field equations and the final line-element 
describing the radiation envelope of a star is 


—_ = —(I- 2m-* r) dr? — r?(d0? + sin? 6dd?) + 7 mt a an) dt?, 
(4-11) 
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i PO Eg Ri MPa) 
x : a er ES st eA 


with m' (1- an) =f(m), m= m(r, t) 
°) for Wy Sr<n, to<t<h. 
es The surviving components of energy tensor are 
m'? m 
: —T=Ti= gg, T= gages T= — gee (4-12) 
6) (Vaidya, 1943). 
5. THE OPERATOR d 
dt 
7 The relation (4-4) 
dm _ 
™ 0 (5-1) 
e 


is a type of relation peculiar to the field we are investigating. In this section 
we obtain some more relations of this type. On eliminating v‘ from v,0" = 0 
n and (v), v0” =0 we find 


dv! dv! a rv’ + vy’ . is 
n &f oF a rer vt( 5 Yt A ew) *) = 0, (5:2) 


But the last term on the left hand side can be shown to vanish by using the 





) field equations (3-8) and (3-10). Hence (5-2) becomes 
, dv} 
a =e (5°3) 


Another such relation can be obtained by starting with the cquation of conti- 
fe nuity (3-4) 


) (p wv), =0 
which when written out in full gives 
ry aint 19 
5x (r? sin Apu" eat > = (). (5-4) 
) &§ When v is eliminated again we find 
D9 (p2ant A—9) 20. Cp2and) 1 Cp2rnt (* ee 24 anna) : 
> pu') +e yp Pe) + pv) ) + AeA-"2 =, (5-5) 
Like (5:2) this also reduces further to - 
, , (rpvt) =0. (5-6) 


(5-3) and ie together can be used to obtain 





_ (r?pv'v1) =O or . (*T™) =0 (5-7) 
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—_— 
and Ps (rp) =0 (5: 8) 
the former of which would again imply T,2=0 as seen in the last section. 
From the definition of the operator 


d _ pt 24 a? 


de ow me? 


it is clear that it differentiates following the lines of flow. Hence the rela- 
tions (5-1), (5-3) and (5-8) show that m, v! and r2p are conserved along the 

- lines of flow. Here we shall try to understand the phrase “‘ conserved along 
a line of flow”. At any time ¢, a spherical wave-front of radius r= the 
radius of the star starts moving onwards. At every point of this wave-front 
the functions m, v', r?p have certain values at the start. The functions will 
retain these values at any point of this wave-front throughout the motion 
of the wave-front. The boundary r=r, at t=t, is a wave-front. At 
t=), this wave-front started moving with radius r=ry. At that time the 
value of m on the wave-front was M the total mass of the star. Our conser- 
vation result now asserts that on this first wave-front, the value of m will 
always be M. Thus at t=4#, on the boundary of the radiation zone r=r, we 
find m=M. 


We may add some simple mathematical properties of the conserved 
results. If 


dt _ 0, 


dr 
in our co-ordinate system, it means 


uy dy —0 
s* axe Dx” 


which means that d#|dx# is a null-vector. 


It follows that 


OY oeg 
(3.6 2 oxe® rane, (5-9) 
For a line-element of the form (3-1), we have from (5-9) that if 
d d 
“= = 0, . (e#’) =0 (5-10) 


which will show that (4-10) is a consequence of (4-4). 
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The actual values of v' and vt may now be deduced. From (5-1) and 
(5:3) we have 
dv? sm’ dv} 


3 me ilies 
Hence 

v'=4(m), = — es 4 (m). (5-12) 
¢(m) is now to be obtained by using any one of the equations 

pipe mM pa Mm oo 
+s Amr?’ Ty 4nmr? ’ MN 4nr? * 

Thus 42r%pv,v! = — m 
or 

{f (m)}? = vod (1 - om) 

4nr*p r 

or 


yi= ey a m' hs 12 } 


\4ar?p) ” m (4nr?p (5-13) 


6. THE BOUNDARY OF THE RADIATION ZONE 


For the field of the radiation zone of a star we have two boundaries: 
(1) the boundary séparating radiation from the material contents (or the 
internal) of the star and (2) the outer expanding boundary of the radiation 
zone separating it from ‘empty’ space beyond. We shall try to find the 
conditions at these boundaries which will ensure a unique solution. 


The line-element under @iscussion is 


-1 2 
PY oe (1 _ a=) dr?— r? (d6? + sin? 0dd2) +s (1 = ma (6-1) 


m (1 _ om) = f(m). (6-2) 


It contains two arbitrary functions. f(m) is one of them. The other is an 
arbitrary function of t, say ¢(f), which appears, when we solve the partial 
differential equation (6-2) for m. 


The expanding bounding surface of the radiation zone has been taken 
to be a sphere of radius r=r, at a time t=+,. Obviously r, and f¢, are inter- 
connected: We shall now say that this bounding surface is a sphere of 
variable radius r= R(t) which would, of course, mean that R(t) =r,. 
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Beyond the bounding sphere of the radiation zone the space is ‘ empty ’ 
and the line-element is 


pall 
dst—=— (1 — MY" ar? — 72 (a0? + sin? dg) + (1 — *™) a, 6-3) 


We show now, that the continuity of g,, at r= R (4) will be sufficient, firstly, 
to locate the boundary at any time ¢, i.e., to determine the function R (a), 
secondly, to find out the arbitrary function 4 (ft) and thirdly to ensure that 
the total energy of the distribution is M. 


Let Vi(m,r) = ¢(t) (6-4) 


be the general solution of the equation 
om 2m 
” (1 - )= f(m) 


the condition for which is 


2m 
ov (1 - , = fem) al (6-5) 
The value of m is given by 
av 
sm = ¢, (6-6) 
Continuity of g,, gives, at r= R (t) 
m=M 
m= — f(M). 
(6:4) and (6-6) then give 
V (M, R)= a (6-7) 
—f(M) Si 2 (6-8) 
dV dV dV aV 
Here yM and oR denote the values of the derivatives os and 5p Tespec- 


tively at m= M, r=R; whichis equivalent to saying that they denote the 
corresponding partial derivatives of V when the variables m and r in V are 
replaced by M and R. (6-7) and (6-8) are the equations to determine the 
two functions R(t) and ¢(¢). 


To eliminate ¢(t) between (6-7) and (6-8) we differentiate (6-7) with 
respect to f, to get 
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which with (6-8) gives 


(6-9) 


We shall now compare (6:9) with (6-5). (6-5) is a relation in m and r. 
When m= M, r=R, it becomes 


aV 2M aV 
RCE oe) =~ SOO Se ls 


Comparing (6-9) and (6-10) we find 


_ 2M, 


R=1-“) 


The general solution of this last differential equation is 
R + 2M log (R — 2M) — t=a contsant, 


which in our former notation, would mean that if the boundary of the radia- 
tion zone is r=r, at a time t=4,, 


r, + 2M log (r, — 2M) — ft, =a constant. (6-11) 


The function ¢(f) is now given by (6-7). It is interesting to note that the 
boundary radius r= R(t) is determined independently of the nature of the 
function V (m, r). 


Before we proceed further let us study the condition m= — {(M) at 
r=R(t). It says that, at all times, on the boundary of the radiation zone, 
m is a constant. But m is not conserved along a line of flow. Using the 
explanations of the last section, we say that the radiating star goes on 
emitting a series of wave-fronts. As m is not a conserved function, it is 
not constant for each one of these wave-fronts. But as m contains an arbit- 
rary function of f, it is possible to select this function in such a way that 
m takes up a constant value on a particular wave-front. And this is what 
we have done by the condition (6-8). Note that the continuity of g,, at 
r=R ensures that at the start, r=ry, t=, the line-eliment is again (6-3). 


The conditions at the boundary r= R(t) have left f(m) undetermined. 
We expect that f(m) will be determined by the conditions at the inner 
boundary of the star. It is clear that f(m) is governed by the conditions in 
the interior of the star, different stellar models giving different forms of f(m), 
That this will be the case, can be very easily seen from the definition of f(m): 


f(m)=m' (1 _ om) 
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or approximately f(m)=m’ or again f(m)= — m, because m’ is almost 
equal to —m. Thus f(m) measures the luminosity of the star, at the 
Newtonian level of approximation. 


Lastly we may now verify that the principle of conservation of energy 
holds good. The line-element (3-1) can be expressed in the form 


a 
ds? == — {(dx)* + (dy)? + (dz)?} — : 73 . (xdx -!- pdy + zdz)? + e” dt? 
(6-12) 
By using the well-known formule? the energy content of (6-12) is found to 


be 
E=lim {$r(e* — 1) e&2} (6:13) 
r>co 


Hence for all distributions for which the line-element (3-1) goes off conti- 
nuously over some boundary to the Schwarzschild’s form (6-3), the principle 
of conservation 
E=-M 
holds good. 
7. THE ELECTRO-MAGNETIC FIELD 


The outside of a radiating star is the seat of electro-magnetic phenomena. 
So the field which we have considered above must be capable of being 
obtained from an electro-magnetic potential K,. That this is the case, 
has been already shown elsewhere (Narlikar and Vaidya, 1947, 1948). 


8. PARTICULAR SOLUTIONS OF THE EQUATION m’ (1 ~< _ = f(m) 
We shall here solve the equation 


om 2m\ __ 
= (1 — >) = f(m) (8-1) 


under different assumptions for f(m). 


Case (1): Let f(m) be a constant. 


f(m)=k < 4. 
m is given by the algebraic equation (m — ar)* (m — Br)? =¢4 (1). 
Here a, B=43{1 + (1 — 8k)"3}; 


A, B=4{1 ¥ (1 — 8k)-¥}. 


¢@(t) is an arbitrary function of f¢. 


1 Formula (91:1) on p. 232 of Relativity Thermodynamics and Cosmology by R. C. Tolman 
(1934) was used, 
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Next let f(m)=k=} 
m is given by the algebraic equation r=4(m — r) {log (4m — r) + 4(H}. 
Finally let f(m)=k > 4. m is given by the equation 


2 tan es '— 1) (8k — 1)-¥* } = (Bk—1)! log (2m*—mr-+kr®)-+¢ (1). 


Case (2): Some more particular solutions may be obtained by trying 
the following method. The condition that the total differential equation 


m'dr +- mdt — dm= (8-2) 
should be exact is given by 
dt _ dr _(r—2m)dm _ (r — 2m) d(m) 
0° 1 4f(m) (8-3) 


7 rm {a + 2f(r — 2my} 


Solving (8-3) we try to obtain m as a function of m, r and ¢. Then this ri 
and m’ from (8-1) will make (8-2) exact. The solution of (8-2) will therefore 
give us the final solution for m. It can now be verified that the following 
is a solution of (8-3) 


m (r — 2m) (m — ar)” (m — Br)” y= 4 (0) (8-4) 
where (a + B)2n=— 1, a”B"y-f(m) =k” (8-5) 
a=a(m), B=B(m), y=y(m), k=a constant, (8-6) 
(eo + >") =k (3 + 2n) (40m) yy" (8-7) 
and f(m) is to be taken as 
f (m) 1 1 +anyt0+” + 2n\'+2” 
(5 + a ) =¢(4f(m) e ie } (8-8) 


n, c are constants andl +n+40,3-+2n0. Various cases follow from 
this solution for different values of n. 


Case (3): 1 +2n=0, f(m)=cm—1, c a constant. Then 
m i i oe an) =cm— 1, 


4 cm? (r — 2m)? = = $2 (cm — 1)? (er? — - 4r + 1 4n). 


The en solution for m is 


2(7-" +2 a * og {(r~ <) + (“-% : 4 +a" 


— 2 Ig (cm — 1)=¢ (2). 
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Case (4): Let 5 + 6n=0 then 50f(m) = (625 cm + 36)"?, ca constant. 


m(r — 2m) ae 7. 9 \ve .— 2 ; 3f 9 \5/6 
f(m) X (m, r)= 4, ( | 200) = X, (m 5 mr [ 8 -{- 200). 


The final (m, r) relation is given by 


- 9V =~, r—2m 
7” f om * a * ) am, 


where the integrand on the right hand side will be a function of m only and 
V= f Xr dr + (ft). 


In the last expression while performing the integration with respect to r, 
m is to be treated as a constant. 


My thanks are due to Professor V. V. Narlikar for having suggested 
this problem and for general guidance during the work. 


SUMMARY 


A star of mass M and radius ry is supposed to start radiating at time fo. 
The zone of radiation extends to r=r, at a later instant t=f,. The energy 
tensor for the radiation zone, describing the directed flow of radiation, is 
evaluated and a relativistic line-element representing the field of radiation 
for r9 < r<r, and corresponding tg <¢ <1, is obtained. It is shown that 
certain quantities m, v1, rp, etc., are conserved in the field along a world-line 
of flow. At r=ro, t=t,) and at r=r,, t=4¢,, the line-element reduces to 
Schwarzschild’s static form for a mass M. The conservation of energy is 


verified. The electro-magnetic potential K, of this field has already been 
obtained elsewhere. 
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Art. 1.—It is well known that Tschebyscheff’s function T,,(z) forms a 
particular solution of the functional equation: 


Sorts (2) ‘ae hn (z) + + u-1 (z) — 0, (n > 1) (D 
as also of the differential equation: 
», d2w dw 
(1 — z?) dz —2— + n?w=0. (A) 


It is interesting to verify that T,,(z) satisfies another functional equation, 
ViZ., 


(1 — 2°) f,'(2) + ntfs (2) — 4 f-1 2} = 0. (n > 1). (IT) 
We shall write (I) and (II) in the equivalent forms, 


bfy-1 (2) = 2h, (2) — : (2° — I) fa (2) (III) 
and fra (2) = AZ) + : (2 — DS, ©). (IV) 


We now propose to investigate all possible common solutions of (1) and 
(II) or what is the same thing (II]) and (IV). We shall henceforth refer to 
the differential equation (A) as the T-equation of rank n and symbolise 
it as T™. 


To fix our ideas, suppose that for some fixed value of n, f,,(z) is an arbi- 
trary solution of the equation T™), so that 


(1—z")-f,"(z) — 2f,'(z) + n°f{z) = 0. (1) 


Then introduce the two “ contiguous” functions ;4,_,(z) and f,,., (z) in accor- 
dance with the relations (III) and (IV). 


Differentiating (I1V) and utilising (1), we find 
urs’ =O + D fa) + pL @}. 
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A second differentiation, combined with (II), leads to: 


2 avr” (2) = * + Ot nit. —_ | I, @)— : ot ‘Sf, 2). 


n n 
(3) 
The relations (2) and (3) accord with the relation: 


( a 2") fier” (z) a —s (z) + (n + 1)? fsa (z) =0, (4) 
showing that f,,., (z) satisfies T+. 


Further, from (IV) and (2) we get: 
l 2 , 
2 [2 fan ()— ‘+ 1 (*— I) fsa | 


=: 2 a) eee nL ae) 
“am y+! —)f, Of- Fp et LO + 


: l 2 , 

1.2., tf, (z) as 2S asi (z) aan n+1 ‘ > ce Dfnsr (z). (3) 
This signifies that f(z) is one of the two functions ‘ contiguous ” to 

f eva (2). 


If we now define the second “ contiguous ” function f,,,. (z) of fn4. (z) by: 
l ‘ ous 
2f t2 (z) as ZS n+i(Z) 2 n + 1 -(z?§— 1) f sett (z), (6) 


and recognise that the other “ contiguous ” function of f,,.,(z), viz., f, (z, 
satisfies a relation, viz., (5), which is similar in form to the equation (III) 
for f,,-1 (z), we infer on the strength of the results proved before that /,,.. (z) 
satisfies T”*». Proceeding in this way, we derive an “ ascending” sequence 
of functions, viz., 


| (Z), Font (z), S wi (z), (7) 
which conform respectively to the sequence of differential equations: 
TM), Teen, T+), 


and between any two contiguous members of which there subsist the two 
functional relations (1) and (II). 


We shall now look for a “‘ descending ” sequence of functions, possess- 
ing similar properties. 


For this purpose we may differentiate (III) and utilise (I), so as to deduce 


tfea O=—- {FO —-f, Of. (8) 
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Another differentiation combined with (1) leads to: 


fia" (z) =a [- = + — ' l =] de (z) _*" n —? Q Ss (z). (9) 


n n 1 — 2 
Manifestly (III), (8) and (9) abide by the relation T'-», viz., 
(J iat 2) Sea (z) een Sa (2) a (n al Fs (z) = 0. (10) 


Repetition of the previous line of argument leads to a “ descending ” 
sequence of functions: 


IAD. $04 OS oa @r: --Se@, (11) 


satisfying the sequence of equations of the type { T’” } and at the same time 
the pair of functional equations (I), (11) or (LD, (IV). 


‘ 


It is scarcely necessary to add that, whereas the “‘ ascending ” sequence 
(7) comprises an infinite number of elements, the ‘* descending” sequence 
(11) comprises only a finite number. 


Art. 2.—It is now easy to make use of the results of the foregoing article 
to investigate the common solutions of the sequence of differential equa- 
tions of the type T” and of the pair of functional equations (I) and (II). 
The needful for this is to take the initially fixed value of n, as contemplated 
in §1, to be unity. That is to say, we have to start with an arbitrary solu- 
tion f, (z) of T, so that 


d—2)-f"O-A@+A@=0. (12) 


Application of the foregoing process at once places at our disposal a 
sequence of functions, beginning with f, (z), viz., 


Si (2); So (2), teeny £ (z), gauatave a . 


which are compatible with (I) and (II) or (III) and (IV). As regards fo (z), 
which is one of the two functions, “‘ contiguous” to f,(z), we have, on 
putting n =1 in (ill), 


So(2) =2{ (1 — 2°): fy (2) + 2h @}- (13) 
In as much as the left side of (12) is 


=7--2) Ot HO} 
we get 

(1 — 2?) - ff’ (z) + 24 (2 =const. (say, c); (14) 
so that 

f(z) =2¢. (15) 
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Thus finally we obtain the common solutions of the pair of functional 
equations (I) and (IL), or (IID) and (IV), viz., 


pa) 4 eS eee 7 See 


We may remark that the common solutions involve two and only two effective 
numerical parameters, these being the two arbitrary constants inherently 
present in the arbitrary solution f,(z) of the differential equation of the 
second order T. It is needless to point out that the constant value (2c) 
of f(z) is not an independent entity, but is rather dependent upon the two 
constants of f,(z) by virtue of (14) and (15). 


It will be observed that the common solutions of (1) and (ID) automatically 
satisfy the sequence of differential equations {T”}*. 


In the succeeding Article we propose to scrutinise in some detail the 
relationship that’ subsists among (III), (IV) and T”. . 


Art. 3.—Starting with the two functional equations (III), (IV) and the 
differential equation T™”, 


we shall combine them in three distinct pairs, viz., 
(111) and (IV), (IIIf) and (A), (IV) and (A) 
and study the legitimate inferences that can be drawn in each case. 


Case I1.—First suppose that f,,(z) satisfies (II1) and (IV), (see Art. 1). 
Then from (IV), 
be Afe + 5 GDL O} Sen’ © 
=F {tin O- 44 O} from (II) 
-hti [et{so+}-@-ps-e} -440)], 


from (IV) 
which when simplified leads to (A). 





* The differential equation T(°) being of the degraded form: 
(1 — 2?) fo”) — zf'o(2@) = 9, 
t is clear that f, (z) = const. is one of its solutions. 
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Case II.—Next suppose that f(z) satisfies (III) and (A). Then 


b{ofn1@ +54 Ser’ O} 
_e \2f, Gi : ‘ir = Df, (2)! 


2 ] “—" P l — Zz a 22 , 
+ — 1 r \tn (2)+S;, (z) + n . tn (z) = n ; I, 

from (LD 
=f,(z), on simplification by use of (A). 


oe Zh n-1 (z) + <4 . f (z) — at (z), 


. ] . oe ; 
or, 2h (z) + nl . (- = lf Nt (Z) — 2h ts (2), 
which is no other than (IV). 
Case IlI.—Now suppose that f, (z) satisfies ([V) and (A). Then 
i 
os : 4 . me 
= 12h wt (2) n+ l a Vy) fia (2)¢ 


l 9 , 
a \fu(2) + n : (2° a li fn ()} 


1 ‘ ° , \ | 1 2 ” 
~ 7 CS -—DIS.L.O+LO +, - VFO 


2z ; : 
v7. Fuh, from (IIL) 
=f, (2), on simplification by aid of (A). 
2 ; 
oe Z Sts (z) aa n | 1 ee = IS ust (z) =f, (2), 


or, 2f ,{2) aD : F (z? a 1) f,,(2) =4tf,,.(2), 


which is the same as (ILI). 


Thus the three equations ({\ly), (1V) and (A) are so related that the 
combination of any two of them automatically ieads to the third. 


Art. 4,—We shall close this topic with a brief reference to certain special 


_ Solutions. It is directly verified that each of the three functions: 


1 l ; 
FA2)= gna» sin (n COS™'z), 51 - cosh (n cosh ~*z), as sinh (” cosh-1z) (16) 
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satisfies the two functional equations (I) and (II) and the differential equation 

(A). Besides, the function V,,(z),* defined by ‘ 
1 1 nn+l1 1 

V,(z — “Fen . 2 , F(5. 3 ’ 1+ ny x) (| z | > 1) (17) 

and called ‘ Tschebyscheff’s Function of the Second Kind” is known to : 

satisfy (1) and (A). So V,(z) musi conform to the equation (II). 


As remarked in the paper referred to in the foot-note, the general solu- 
tion of (1), which is after all a functional equation of the second order, can 
be thrown into the form:+ 

f A2) = 4,(2)°8 (2) + By (2)-$ (2), (18) : 
where a, (z) and £,,(z) are two particular solutions of (I) and @(z), 4(z) are 
two arbitrary functions. 


For obvious reasons, a, (z) and 8, (z) can at pleasure be replaced by any 
two of the functions T,,{z). V.,(z) and those marked (16). 











* Refer to Bagchi and Chakravarty’s paper, viz., ‘‘ Note on Tschebyscheff function T,. (z) and 
its associated equations,” where certain properties of the ‘‘so-called’’ Tschebyscheff’s function 
of the second kind are discussed. (Vide Journal of the Indian Math. Society, 1950.) 
+ Throughout this paper, the parameter n is supposed to be a positive integer. 
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1. INTRODUCTION 


BARITE or heavy spar is a well-known minera! of which colourless and trans- 
parent big-sized crystals are easily obtainable. Hence it is not surprising 
that it has been a subiect of investigation by various workers in different 
fields of study. Thus, R. S. Krishnan (1946) has studied its Raman spectrum 
with some thoroughness; while the temperature dependence of Raman 
frequencies and their width has been investigated by Narayanaswamy (1948), 
Amongst the directional properties, mention may be made of the elastic 
constant determinations of Voigt (1910) by static method and of Bergmann 
(1938) by dynamic method. Bridgman (1928) has measured the linear 
compressibilities of this crystal along the crystallographic axes. Results 
of magnetic susceptibility (K. S. Krishnan and others, 1933) and dielectric 
constant measurements (Rao, 1949) have also been reported in the literature, 
More recently the thermo-optic behaviour of this crystal has been studied 
in detail by Radhakrishnan (1951) in this laboratory. 


In view of such diverse studies of the properties of this mineral, it is 
surprising to note that its thermal dilatation has not been studied so far. In 
LC.T. (Vol. III, p. 44, 1928) is reported one value for the coefficient of 
linear expansion which appears to be misquoted. As an accurate know- 
ledge of the directional dependence of thermal expansion is required for an 
understanding of the relation between thermal behaviour and crystal struc- 
ture, the author has undertaken the measurements of thermal expansion 
of this crystal in different directions and the results of study are set forth in 
the present paper. 


2. PROCEDURE OF THE EXPERIMENT 
Barite is an ionic crystal belonging to the .bipyramidal class of the 


orthorhombic system. The relevant crystallographic data for identifying 
the faces are given by Groth (1908) according to whom the axial ratios are 
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a:b: c=0-8152:1:1-3136. The c-face (001) is a perfect cleavage plane. 
Cleavage along b (010) is also perfect, while along m (110) it is rather good. 
The symmetry of the crystal under investigation being orthorhombic, the 
three axes of the thermal ellipsoid coincide in direction with the crystallo- 
graphic axes, and thus there are only three independent coefficients of 
expansion, viZ., a4, G2 and a33. Hence, for the purposes of the present 
investigation, three sets of spacers, in the form of right pyramids, with their 
altitudes coincident with the crystallographic axes, are needed. Two more 
sets of pyramids with their altitudes along (i) a direction equally inclined 
to the three axes and (ii) along a direction perpendicular to the natural 
cleavage m-face have been employed for verification and a check on the 
values previously determined. 


The thermal expansion coefficient in an arbitrary direction is easily 
obtained by a transformation of axes to the concerned direction. The new 
set of axes X,’, X,’, X;’ are defiued with respect to the old ones by the usual 
direction—cosine scheme :— 























| XY) Xe | | 

| | 
Xi | Ci Ci2 | Cis 

| 
| 

Xe | Ca C22 | Coz | 

| 

Me Ca | Css | Ces | 





As usual a’;, indicates the constant in the arbitrary direction, which 
coincides with the X;’-axis. The expression for a3,’ is given by (Wooster, 
1936)— 


33° = Cy37a11 + C4370.99-+ 2Cj3Cog12+ 2C3C33413-+ 2Co3C ages (1) 


in the orthorhombic system a; (ij) vanish. The direction cosines for 
the arbitrary direction (i) are given by 


1 
Cu = Coo= C33 = V3 
O33 = $2 a; js (2) 


For the expansion perpendicular to the m-face, the expression (1) further 
simplifies 


Hence we get 


Qian, = Cy32a,)+ Cy37ag0. (3) 
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The specimen of the crystal of barite at the disposal of the author had 
its cleavage planes c and m perfectly well marked out. The axes were 
located with the help of ordinary polaroids. The c-axis was contained in 
the m-plane, while ‘a’ and ‘b’ in the c-plane. The latter were identified 
by the length of intercepts made by the m-plane with them. 


3. RESULTS AND THEIR DISCUSSION 


The coefficients of thermal expansion, a,, a3, a, along the three crystallo- 
graphic axes a, b and c are entered in Tables I-III. The calculated values 
of the same from the interpolation formule are also put forth in the same 
tables. 

TABLE | 
Coefficient of linear expansion for barite parallel to the a-axis 
ag = 0-0,1362 + 0-0,1298 ¢ + 0-0,91192 














Temperature | aq xX 10° aq X 10° | Difference 
“Cc | observed calculated pA | 
| 

44-3 14-23 14-22 2 

78°5 14-67 14-71 - 3 | 
111-2 15-24 15-21 +2 | 
142.5 15-83 15-71 +-7 | 
175+6 16-26 16-26 “ | 
210-6 16-76 16-88 - +7 
244-3 17-48 17-50 = 
276+4 18-27 18-12 +8 | 
307-3 18-83 18-73 + +5 
340-8. 19-36 | 19-43 - 4 | 
377-1 20-21 20-21 ; | 

TABLE II 


Coefficient of linear expansion for barite parallel to the b-axis 
ay = 0:042395 + 0-0,1256 tf + 0°0491368 7 














| | 
Temperature as X 10° ag x 10° ‘ 
"C observed calculated Difference 
42-6 24-58 24-51 | + -3 
74-5 } 24-98 24-96 oe 
105-6 | 25-48 | 25-43 xi 
136-1 | 25-86 | 25-91 | =—« 
170-9 | 26-42 26-50 | = 
209-6 | 27-18 27-10 | oe 
F 247-0 | 28-01 27-89 + -4 
283-5 28-68 28-61 | + -2 
319-1 29-29 29-35 | we ff 
353-8 30-12 30-08 | oo 
30-87 8 
| 








286 S. S. Sharma 


TABLE III 


Coefficient of linear expansion for barite parallel to the c-axis 
ae = 0-041407 + 0-0,1520 t + 0-0491102 z* 

















Temperature | ac X 10° ac X 10° ete 

ba | observed calculated Difference 
45-5 14-73 14-78 Low ly 

| 10-8 15-09 15+20 | wm 

| 95-2 15-62 15-62 as 
118-7 16-08 16-03 [| 
141-4 16-51 16-44 | fw 
167-1 17-00 16-92 | 4 4g 
195-5 17-44 17-46 | a 

223-2 18-03 18-01 : 
249-8 18-67 18+55 + +7 
275-6 19-09 19-09 ‘i 
303+7 19-72 19-70 A | 
334-1 20-23 19-38 ae | 
363-6 20-97 21-05 as Te 

| 





The values of a,,,, the coefficient of expansion along a direction 
equally inclined to the axes are entered in Table IV. The interpolation 
formula 

@,5- =0-0,1720 + 0-071369 ¢t + 0-0,91204 2? 
has been derived from the experimental values of a,, a; and a, by employing 
the relation (2) given above. The calculated values of a,,. from this 
formula are entered in the third column of this table. From a perusal of 
these values it can be seen that the agreement between the experimental and 
calculated values is quite satisfactory. 


TABLE IV 


Coefficient of expansion for barite in a direction equally 
inclined to the axes 








Temperature Gabe X 106 Gak X 10° 
*e. (observed) =4 3 a;; 
! eee! See ee 

48+4 17-91 | 17-89 
719-2 18-31 18-36 
109-0 18-81 18-84 
138-0 19-28 19-31 
166-1 19-77 19-81 
196-7 20-36 20-36 
229-7 20-94 | 20-98 
261-6 21-63 21-61 
292-6 22-24 22-24 
322-6 22-89 j 22°87 
351-9 23-41 23-51 
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The values of a,,, are entered in Table V. 


TABLE V 


Coefficient of expansion for barite perpendicular to m-face 
a = 0°041770 + 00,1254 t + 6-0491262 2? 



































| | 
Temperature | a, x 10° a, X10° | Difference 
C. | observed | calculatec| | % 
~ oy ~ sae 
43-4 18-28 18-27 Re 
} 75-1 18-84 18-71 + -7 
105-7 19-28 | 19-17 + -6 
| 135 +5 19-72 19-63 + +5 | 
164-6 20-07 20-10 —- 2 | 
195-8 20-59 20-64 - 3 | 
228-8 21-23 | 21-23 a 
260-8 21-94 | 21-83 + +5 
291-8 22-55 | 22-43 +5 | 
321-9 23-04 23-04 . | 
355-0 23°62 | 23-74 — <4 | 
' 
32, 
y 
> 
} 
7 
| 28} oh 
2a 
‘oO. 
2 
y 20h 
16F 
12 , + 
oO 80 160 240 329 FOO 


Tempera Eure C.—+> 


Fic. 1. Variation of coefficients of expansion of barite with temperature. 
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The variation of the various coefficients of expansion with temperature 
has been depicted by smooth curves in Fig. 1. 


a,,, can also be evaluated from the observed values of a, and a, by 
employing relation (3) given above. Thus if @ be the angle between the 
normal to the m-plane and a-axis, the relation (3) reduces to 


a=a, cos? 6+ a, sin? 0 (4) 
where 6 = tan“ 0-8152 = 39° 9’. 


For the sake of comparison the experimental values of a.,, are tabulated 
below along with those calculated from (4). 

















TABLE VI 
Coefficient of expansion (Xx 10°) perpendicular to m-face 
| 
Temperatute °C, oc. | 100°C. | 200° C. | 300° C. 400° C. 
| 
om | : ) 
a,,, (observed) 17-70 19-08 | 20-71 22-60 | 24-74 
| | s cscahaeapl 
a,,, (calculated) 17-74 | 19-14 20-80 22-72 24-89 
Equation (4) | 























Although there is a systematic deviation between these two sets of values, 
yet the discrepancy is not very marked and the agreement can be taken to 
be fairly satisfactory. 


The values of the principal coefficients of expansion at 0° C. have been 
calculated and are given below:— 


a, = 13-62 x 10-*; a;=23-95 x 10-*; a,=14-07 x 10-6, 


It is seen that the, thermal anisotropy is maximum along the b-axis, while 
along the a- and c-axes the coefficients are nearly equal. The exact in- 
equality can be represented by a, >a,>a,. It may be mentioned here 
that in this respect the thermal behaviour of barite is different from its optical 
behaviour. According to Groth (1908) the principal refractive indices are 
n, = 1-6491, ns=1-6381 and n,=1-6369. Thus in this case n, >; > n,. 
The optical and thermal ellipsoids are differently oriented in the crystal. 
But a survey of the literature shows that many other properties of barite, 
namely dielectric constant, magnetic susceptibility, linear compressibility, 
etc., exhibit anisotropy in the same sense as the coefficient of thermal 
expansion, i.e., being maximum along the b-axis and minimum along the 
a-axis, 
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In conclusion, the author expresses his heart-felt thanks to Professor 
R. S. Krishnan for his constant interest in the work. 


SUMMARY 


The paper describes the results of study of the thermal expansion of 
barite and its directional dependence. The principal coefficients of linear 
expansion along the crystallographic axes are given by 


a, = 0-0,1362 + 0-0,1298 ¢ + 0-0,91192 #2 
a; = 00,2395 + 0-071256 ¢ + 0-0,)1368 7? 
a, == 00,1407 + 0-071520 ¢ + 0-0,91102 2°. 
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INTRODUCTION 


THE object of the present note is to trace the connecting link between the 
Tschebyscheff’s function (“of the first kind”), viz., T(z) and the 
Tschebyscheff’s function “of the second kind”, viz., V,,(z) introduced in 
a paper by H. Bagchi and N. K. Chakravarty, published elsewhere.? 


1. Defining the sequence of constants {c,,,} by the relation: 


1 


_ (Tv) do 
Co, p= sa v2)t ’ 


—1 


where T, (v) is Tschebyscheff’s function of the first kind and the parameters 
p, q ate positive integers, we note the relations?: 
(i) ¢,, s=9, if q<p, 
7 


and (ii) cy »= 2-1 if g=p. 


By aid of a suitable formula of reduction, one can easily substantiate 
the following additional properties : 
(iii) c,,,=0, if gq > p and q — p be an odd integer; 
and (iv) Ce. p= san v*’C,, if q>p and q—p be an even integer 2r. 
In particular, we note the following special cases of (iv):— 
e a#+2. a  (n+3)(n+ 4), 


Cyt2, 2 = F2n-1 © 7° Cnt4, n= 22n+3 2! ‘ 


a  (n+4)(n 3 5) (m+ 6) and so on. 


Cut, n 2275 


290 

















Mutual Relation between Two Kinds of Tschebyscheff’s Functions 291 


2. Tschebyscheff’s function V,,(z) ‘“‘of the second kind”, has been 
defined by 


Vv an,! ; E . ree .f 4 n(n+ I(n+2)(n+3) 1 
” 2” (2%  41!@4+1) 7? 4-2'(n+1)(n+2) = 7A 
4 aint Yat+2)(n+ 3)\(n+ 4) (n+ 5) | 1 \ 
43-3! (n + 1) (n + 2)(n 4+ 3) pt oo] (1) 
where |z|>1. 
Plainly the R.S. 
=p - [4 gat l eels, eo ] a 
2m Lne® 22° (nt+22 2! 71 a+ eae Oe, 
_(n+3)(1+ 4) 1 | __, (n+4) (n+5) (n+6) 
=z (n +6)z"+6 «22745 Z3 "h +e+-to = | 


n 1 1 1 
on : [a Cron + (n + 2) z*+2 Cote, n + (n + 4) z7t4 * Cota, a 


1 


+ (n +6) ztt6 * Coté, 2 Pass } by Art. 1. 


Using the relation: 
Co,n=9, if g #n and q—n be odd, 


we can fill up the missing gaps in the sequence of terms in the above series 
and thus write 
1 ] 1 ] 


n 
V,{z) rs an : [. "Ce t 272 Cent 373° nm beree $ (n —1)z = ee 


1 1 ] 
+ nz” Chm + (n > 1) 2") Cath, n tee 2) ze Cr t2, at +++. to oo | 


1 


nr @ fi (v)'T,(v) 
“ow Ja oe pte. (2) 
—-1 





Inasmuch as the series on the R.S. of the familiar equality: 


—log(1 —*)="F 1 - (2), (—1<v<l, |[z|>1) (3) 


r=1 
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is absolutely and uniformly convergent for unrestricted variation of v in the 
domain (— 1, 1) and the multiplication by the bounded quantity,* viz., 
T, (2) 
(i — 0% @ 


does not affect the uniform convergence, we conclude on the strength of 
known lemmas of Analysis that the series (3), when multiplied by (4), must 
retain the property of uniform and absolute convergence so as to admit of 
term-by-term integration. 


Consequently we obtain: 


ef ve a (2 (z%- v7) de me - Sat OO log, (i _ *) dv. 


Hence by (2), 


V,(@)=— af Gsm toe. (1-2)d, (21>) 


An equivalent formula is: 
n L* 
V.@=—3 Ja nO log (2 —v) do, (|z|>1) © 


(5) or (6) may be regarded as the analogue of Neumann’s integral formula 
for Q,,(z) in terms of the Legendre function of the first kind P,, (z). 


3. The legitimacy of differentiation (w.r.t. the parameter z) under the 


integral sign in the R.S. of (6) of the foregoing article being kept in view, 
we immediately deduce: 


‘io ee ee 
Va (= — 27 fi — 2) (1 ~ oy (1) 
Writing T,,(v) in the form: 


Re (v) = Ba (z) os {T,, (z) am 1 (v)} 


* Strictly speaking, the function (3) is bounded in the interval (— 1+ e, 1 — e), where 
e>0. At either end-point (— 1 or 1), (e = 0), the function (3) becomes infinite and the definite 
integral (4) becomes improper. It can be proved by-a’ priori reasoning that this improper integra, 


is uniformly convergent at each of the two end-points, so that the term-by-term integration is 
valid, 
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and substituting in (1), we get: 


’ n i d: 
V, (z) = Vn Be (z) f (z —%) i —y2)t a K, (2) 
where 
ai F * (z) oe T,, (v) 4 
K= J giao @) 


Since T,, (v) is a polynomial of the type: 
T,, (v) = ayu”+ agv”-? + ag**+ ...... y 
the quantity 
T,, (Z) — Ty (2) 


z—v 
is a polynomial of order (nm — 1) in z, viz., 
a, (v) + 2-1 + ag (v) 27? + ag (v) 27-3 + ...... . 


the coefficients of the different powers of z being polynomials of order not 
greater than n — | in v. 





Hence (3) may be re-written as: 


“K=b,2"1 + bz? + byz"-3 + 


a, (v) ; 
a? dv. 


Now the definite integral, occurring in (2), viz., 


where b,= 


; dv 
(2—v) (1 =v 


when evaluated by the substitution w=(j F $f is easily seen to be 


SRST a ee 


ru 


4 (z? ae 1)t ° 
Hence the equation (2) may be presented in the form: 
oe 
¥e (z) hal 2 : (z? come i R,4 (z), (| Z | > 1) (5) 


where R,_, (z) is a polynomial of order (n — 1) in z, 
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The relation (5) may be looked upon as the analogue of the corres- 
ponding relation between the two kinds of Legendre functions, viz., 


Q,(z)=4P, (2)log 2+! — w,_, (0, 


z—1 
where W,,_, (z)® is a polynomial of order n — 1. 


ArT. 4. We shall now close this topic with a short note on the analytic 
continuation of the function V,, (z). 


For obvious reasons the function V,,(z), as defined by the negative 
power-series (1) of Art. 2, exists in the region, exterior to the unit-circle 2, 
viz., |z|=1. But the domain of analytic existence of the function V,, (z) 
as defined by (6) of Art. 2, is comparatively wider. For the line-integral 


“day (1) 
= 

(taken along the segment A of the real axis, intercepted between the end- 
points 1 and — 1) is easily seen to represent a function of the parameter z 
(real or complex), analytic anywhere in the finite part of the plane except 
at points on A. In fact, the finite segment A (— 1, 1) is the “ so-called ” 
Hermite’s Section for the line-integral (1) and is accordingly a line of 
discontinuity for the function of z, as defined by it. It follows automatically 
that the same rectilinear segment A counts as a line of discontinuity for the 
function V,,(z), as defined by (6) of Art. 2. In other words, the function 
V,,(z), aS defined by (6) of Art. 2, exists at allinterior points (of 2), other 
than those situated on the diameter A. Thus the passage from (1) of Art. 2 


to (6) of Art. 2 implicitly involves the process of analytic continuation of 
the function V,,(z). 
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SYNERESIS OF SODIUM OLEATE GELS IN 
ORGANIC SOLVENTS 


Part I. Effect of area of gel-surface on the syneresis of 
sodium oleate gels in pinene 


By MATA PRASAD, F.A.Sc. AND V. SUNDARAM 


(From the Chemical Laboratories, The Institute of Science, Bombay) 


Received July 20, 1950 


HoLMEs and co-workers! have found that the syneresis of silica gels varies 
directly as the area of gel-surface exposed for syneresis. However, Ferguson 
and Applebey? find that syneresis in silica gels is independent of the surface 
exposed. Chapman and Buchanan* have observed that the syneresis of 
starch gels prepared in larger tubes is greater than that of the gels prepared 
in smaller tubes. It is evident that the study of the effect of surface on the 
syneresis of gels is important and has not so far received much attention. 
The present investigation deals with the study of this effect in the case of 
gels of sodium oleate in pinene. 


EXPERIMENTAL TECHNIQUE 


The following series of experiments were carrried out to study the effect 
of surface area on the syneresis of sodium oleate gels in pinene. 

Series A.—Three flat-bottomed dishes of diameters of (1) 4:3 cm., 
(2) 7-1 cm., (3) 9-8 cm., respectively, were used. It was so arranged that 
the amount of the gel in each dish was proportional to its surface area. The 
thickness of ihe gel layer was practically kept constant. 

Series B.—Varying amounts of the gel were taken in the same dish of 
4-3cm. diameter. 

Series C.—The same amount of the gel-forming solution was taken in 
each of the three dishes used in Series A. 

For all the experiments a solution of 1% sodium oleate in pinene was 
prepared by dissolving the requisite amount of sodium oleate in pinene at 
140° C., according to the method followed by Prasad and co-workers.‘ 

For the experiments in Series A, approximately 5c.c. of the solution 
were transferred to dish (1), 14c.c. to dish (2) and 27¢c.c. to dish (3) and 
then allowed to set in the dishes. For the experiments in Series B, approxi- 
mately 5c.c., 10c.c., 15c.c., 21 c.c., and 27 c.c. of the solution were transferred 
to dish (1) and allowed to set. For the experiments in series C, approxi- 
mately. 27c.c. of the solution were transferred to each of the three dishes 
and allowed to set. All the dishes were kept in a thermostat maintained 
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at 30° +0-1°C. In all the cases the amount of syneresis was measured 
three hours after the gels had set by the method followed ty Prasad, Hatti- 
angdi and Mathur. 
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RESULTS 


The results obtained are given in Table I in which the following nota- 
tions are used: 


D_ diameter of the dish in cm., 


h thickness of the gel layer in cm., 

er ern eee mD (h + D/2) 
Ss Specific surface in sq. cm. per c.c. [ equal to ~~ aD2H/4 |. 
W weight of the gel in grams, 


x amount of syneresis in grams, 



































Ys percentage syneresis equal to ed 
TABLE I * 
Series | D h P w | x Y y’ 
l 
| 453 0-3521 6-609 | 4-356 1-196 | 27-46 
A | 71 0.3407 6-433 11-500 | 2-806 | 24-35 
| 9-8 0 +3460 6-187 22. 250 5-400 | 24-56 . 
| 453 0+3521 6-609 | 5-356 1-196 97-46 | 27-46 
4:3 0-6869 | 3-842 8-500 1-912 | 22-50 15-97 
B | 4:3 | 0-9924 | 2-946 | 13-284 | 2-088 | 17-00 | 12-25 
| 43 | 21-3980 | 2-360 | 17-306 | 2-150 12-43 | (9-81 
| 43 1-8520 2-011 | 22-920 | 2-481 | 10-82 | 8-36 
| 9-8 | 0-3460 6+187 22.95 | 5-400 | 24-56 | 24-56 
Cc | 7-1 0-6651 3-568 | 22-45 4-100 | 18-26 14-18 
| 4-3, | 18520 2-011 22-92 2-481 | 10-82 | 7-99 








DISCUSSION OF RESULTS 


It is seen from the above table that (i) if the specific surface is constant, 
the percentage syneresis is also constant, provided, however, the thickness 
of the gel is small and is of the order of 0-35cm. (cf, Series A), (ii) the 
percentage svneresis depends upon the specific surface, being greater, 
greater the specific surface (cf., Series B and C), and (iii) the amount of 
syneresis per unit surface depends upon the thickness of the gel, if it is small, 
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d and is nearly independent of the thickness if it is large. The values of per- 
i- centage syneresis (Y'), given in the last column of Table I, calculated on the 
. basis of the specific surface, using the value of Y for the gel of the smallest 
thickness in the B and C Series, show that there is no direct proportionality 
- between s and ¥. This observation is in contradiction to the conclusions 
of Holmes and co-workers! which are not supported even by their own 
experimental results. It is clear that the syneresis of gels of fair thicknesses 
is not wholly determined by the specific surface; presumably there are other 
factors as well which influence it. 


The effect of surface area exposed on the syneresis of these gels can be 
explained as follows: 


The number and the size of micelles in a gel depend upon the 
concentration of the dispersed phase and the temperature to which the gel- 
forming solutions is allowed to cool. Now, at any temperature, if a fixed 
volume of a gel sets in vessels of increasing surface area, the number of the 
a micelles at the surface and therefore the number of interfibrillary spaces on the 
surface of the gel will increase. Since the exudation of the syneretic liquid must 
‘ take place ultimately through the interfibrillary spaces on the surface of the 
gel, the increase in the number of such spaces would facilitate the flow of the 
exuded liquid and hence the amount of the syneresis in the same gel in a given 
time would depend upon its surface. However, if the height of the gel is 
increased within small limits an increase in the amount of syneresis may be 
expected, on account of (i) the increased hydrostatic pressure exerted by the 
increased column of the gel and (ii) the possible decrease in the stability of 
the gel on increase in the thickness, provided the increased height does not alter 
the gel-structure or change the number and size of the fibrils present therein. 


SUMMARY 


The effect of the area of gel surface exposed on the syneresis of sodium 
oleate gels in pinene has been determined. It is observed that syneresis is 
directly proportional to the specific surface per unit volume of the gel, pro- 
‘ vided the thickness of the gel is small, and nearly constant. When greater 
: thicknesses of the gel are considered syneresis is slightly higher than the 
value calculated on the basis of the specific surface. 
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INTRODUCTION 


BREYER AND GUTMAN! have evolved a new and interesting technique for 
studying the electrode kinetics of the reversible depolarising processes taking 
place at a dropping mercury electrode by superimposing an alternating field. 
Systems involving the discharge of cadmium, thallium and zinc ions have 
been studied by them with an alternating field of 50 cycles per second. The 
discharge of lead ions has been studied by a similar technique by Doss 
and Agarwal? using the oscillograph for measuring the alternating current. 
They have extended the study to different frequencies of the field ranging 
from 12-5 to 600 cycles per second. The present paper reports the work 
on the discharge of nickel ions studied by a modified technique. 


EXPERIMENTAL 


The circuit diagram is given in Fig. 1. The source of alternating current 
is a B.S.R. oscillator capable of giving frequencies ranging from 0-16,000 
cycles per second. 10 Volts from the oscillator is fed on to a potentiometric 
arrangement ACDB consisting of the resistances R,, R, and R3. This 
arrangement serves to give 45 mv. across R, which is superimposed on the 
dropping electrode along with a D.C. voltage obtained from a Cambridge 
pH potentiometer (with the galvanometer shorted off) used as a potential 
divider. The voltage output of the potentiometer across YZ can be varied 
at will, the resistance R, being 225 ohms per volt of D.C. voltage. 


When the pulsating (A.C. superimposed on D.C.) potential is incident 
on the dropping electrode, a pulsating current is produced. Our present 
interest is to measure the A.C. component of the pulsating current. For 
measuring this, the arrangement used is as follows :—A valve amplifier system 
(Fig. 1) working on the non-linear portion of the characteristic curves of 
the valves, partially rectifies and amplifies any incident alternating potential. 
This amplified current is incident on a Leeds and Northrup galvanometer 
(period: 14-5 sec.; resistance: 1,182 ohms; shunt: 99 ohms). The tappiag 
key K connects the amplifier system to the potential drop across R; in the 
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rest position P, and to that across R; when the key is depressed to the posi- 
tion Pz. First, the resistance R, is reduced to zero, thereby making the 
input to the amplifier system zero. The current passing through the galvanc- 
meter under these conditions is compensated by the potentiometer. Then 
the resistance R; is increased to 3 ohms, thereby making the input to the 

































































Fic. 1 

R = 6002 . £ = 0:25yF 
R, = 13,000 2 C, =O1sF 
R, = 582 LT. =1-5V 
R, =0-3.Q2 H.T. = 45V 
R, = 225 Q per volt of D.C. S =992 
R; = 0—2,500 {2 (according to Re =5MQ 

the current) R, =1MQ 


amplifier equal to 2:3mv. The galvanometer gets deflected. The reading 
on the galvanometer scale is noted. Any current passing through the drop- 
ping electrode passes through the resistance R;. The potential drop thus 
obtained across R; is applied to the amplifier system by depressing the tapp- 
ing key to the position P,. In this position R; is adjusted so as to give the 
same deflection of the galvanometer as was obtained for a potential of 2-3 mv. 
(i.e., in the position P, of the tapping key). By knowing the value of the 
resistance R;, the current can be calculated. The, condenser serves to filter 
off the D.C. component of the pulsating current. 


This arrangement for measuring the current has the following 
advantages: (1) The measurement is not affected by any changes with time 
in the characteristics of the amplifying system. (2) It is also unaffected by 
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the variation of the characteristics of the amplifier with the change in frequency 
of the A.C. 


The dropping electrode has the following characteristics: 
m: 0-69 mg./sec. 
t :5-3sec. in 0-1 M — KCl,,. in open circuit. 


The experiment is first performed with IM KCI containing no other 
dischargeable ions. Hydrogen, obtained from a Kipp after purification 
by passing through mercuric chloride, lead acetate and distilled water is 
passed in the solution for 15 minutes to remove traces of oxygen. The 
A.C. passing through the dropping electrode at the different D.C. potentials 
are measured as explained above and the results are represented in Fig. 2. 
The experiments are carried out at different frequencies of the A.C. 


A.C. curRENT in L(A 








° 0-4 


D.C. POTENTIAL (in vous) 
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The experiment is repeated with IM KCI containing 3-27 x 10-°M 
nickel ions (in the form of nickel chloride) and the results are represented 
in Fig. 3. 
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CURRENT IN ACA, 





A.C. 





O.C.POTENTIAL (iwvoits) 
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Fig. 4 represents the polarographic curve of nickel obtained with the 
same dropping electrode. 


DISCUSSION 


Our results confirm in general the observation made by Breyer and 
Gutman! that a maximum in the alternating current occurs nearly at the 
half-wave potential of the system in question. Fig. 4 shows the polaro- 
graphic curve from which the half-wave. potential comes out to be at about 
1-1 volt. The alternating current maximum from Fig. 3, comes out to 
1-2 volt, which is nearly the half-wave potential. The D.C. potential at 
which the a.c. maximum occurs is independent of the frequency of the a.c. 
field. This is in agreement with the observation made by Doss and Agarwal 
(loc. cit.). There is one important difference, however, between these curves 
and the curve obtained by Doss and Agarwal. Whereas the value for the 
alternating current at the maximum in their system was nearly independent 
of the frequency, with nickel it is found that there is a large stepping up of 
the a.c. current as the frequency is changed from 50 to 100 cycles per second. 
This difference in behaviour can be explained as follows: 
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CURRENT IN A. 





Oc. 





4 og 16 


9.6. POTENTIAL (in voits). 





Fic. 4. Polarographic current : Nickel 


The dropping mercury electrode in contact with the solution of the 
electrolyte presents an electric double layer which behaves like a capacity 
in absence of any reducing ions. The increase in the a.c. current with 
increase in frequency (Fig. 2) is due to the diminution of the capacitative 
impedance. The effect of the presence of the reducible ions is equivalent 
to shorting the electric double layer capacity by a resistance, or more 
generally by an impedance consisting of a resistance and a capacity in 
parallel the value of which decreases with the increase in the intensity of the 
discharge of the reducible ions. The effect of this is to reduce the overall 
capacitative impedance. At the half-wave potential of the reducible ions 
the discharge is maximum and hence there is a maximum in the alternating 
current. It appears that in the system dealt with by Doss and Agarwal, 
this maximum capacitative impedance becomes a very small fraction of the 
total impedance of the system so that the change of frequency has no 
marked effect on the overall impedance of the system at the a.c. maximum. 
With the system dealt with in this paper, however, the minimum capacitative 
impedance is yet comparable with the total impedance of the system and 
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hence an increase in the frequency has a marked effect on the value of the 
a.c. current at the half-wave potential. 


Breyer and Gutman have worked out a quantitative theory to explain 
the form of the curves. One of the assumptions that they have started is 
that the alternating current at any instant would be equal to the polarogra- 
phic current corresponding to the potential incident on the electrode at that 
instant. We have worked out the theoretical curve relating the a.c. current 
with the d.c. voltage on the basis of the above assumption alone. This is 
shown in Fig. 5 (Curve I). The experimentally observed curve (Curve II) 
is given in the same figure for comparison. The above assumption of Breyer 
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and Gutman approaches nearest to the experiment at near the half-wave 
potential for this system at the frequencies of 12-5 and 25 cycles per second. 
At other frequencies and potentials, however, there are large deviations, 
the experimental values being much higher than the theoretical values. This 
divergence is to be traced to the fact that the measured polarographic 
current is an average figure whereas the a.c, current is caused by the 
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comparatively instantaneous changes in the current when the potentials are 
altered quickly. 


A quantitative treatment of the observed phenomena appears to be 
possible only after more detailed and extensive experimental data on such 
systems become available. 


The authors wish to thank Prof. J. M. Saha, Dirtctor, Indian Institute 
of Sugar Technology, Kanpur, for his kind interest in the work. The authors 


are grateful to the Uttar Pradesh Scientific Research Committee for a grant 
which made this work possible. 


SUMMARY 


A method has been described for measuring small alternating currents. 
This has been applied for studying the a.c. currents induced in a dropping 
mercury electrode system containing nickel ions when subjected to a pulsat- 
ing field at differer.t frequencies of the field. This system has revealed interest- 


ing features which are different from the systems investigated by the previous 
workers. 
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THE effect of time on the syneresis of silica gels seems to have been first 
studied by Okatov! who suggested that syneresis is a good measure of the 
maturity of these gels. Ferguson and Applebey* have observed that in the 
initial stages, the velocity of syneresis of silica gel follows an autocatalytic 
curve. Prasad, Hattiangdi and Mathur® have studied the kinetics of the 
syneresis of sodium oleate gels in pinene and have found that, for the interval 
of time studied, syneresis follows a relation similar to the imbibition of 
liquids. The present paper deals with the investigation of syneresis of 
sodium oleate gels in pinene and in xylene carried out over an extended 
period of seven days. 
EXPERIMENTAL 


It has been pointed out in the previous paper* that syneresis varies with 
the surface of the gel. Hence in these investigations test-tubes of the same 
internal diameter, 1-70cm., were used throughout. Gels containing 
0-10 g. of sodium oleate in 10 c.c. of each solvent were used for the measure- 
ment of syneresis. The gels in pinene were prepared in the manner de- 
scribed in Part I. The gels in xylene were prepared in a similar manner by 
preparing the solution of sodium oleate in xylene at 130°C. All the gels 
were kept in a thermostat maintained at 30°+0-1°C. The amount. of 
syneresis was measured by the method of Prasad and co-workers® at different 
known intervals of time after the gels had just set. 


RESULTS 


The results obtained are given in the following table in which the 
following notations have been used: 


t =time after the gel had set, and 
X =the amount of syneresis in grams, 
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TABLE I 
Gels in pinene Gels in xylene 
Zin hrs. 
X | Km X 108 p3 Km x 10° 
1 0-600 20-250 | 1-703 | 60-900 
2 0-795 13-550 1-798 32-370 
3 0-923 10-630 1-844 22-200 
5 1-088 | 7-577 1-861 13-460 
7 1-165 5-823 1-872 9-658 
18 1-317 2-587 | pe Bs 
24 1+374 2-031 | 1-945 2-941 
48 1-485 1-107 2-008 1-629 
144 aa | ea 2-074 0-529 
168 1-729 | 0°375 | 2-090 0-457 








DISCUSSION OF RESULTS 


It is seen from Table I that the gels in xylene synerise more rapidly than 
those in pinene. The graphs of the amount of syneresis against time are 
smooth rising curves (Fig. 1) and are not S-shaped as observed by Ferguson 
and Applebey.* This shows that the process of syneresis of sodium oleate 
gels both in pinene and in xylene is not autocatalytic. 
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Lipatov® has observed that the velocity of syneresis in geranin gels 
follows the law of first order reaction, Prasad and co-workers® have found 
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that the unimolecular constants (K,,) calculated for different intervals of 
time in a given series are not constant but decrease as the time interval is 
increased, thereby showing that the law of unimolecular reaction is not obeyed 
in this case. They have also noticed that the graphs obtained by plotting X/t 
against K,, are straight lines, showing that the relation X/t = KK,, + C, 
where K and C are constants, is fullowed in this case. Similar results have 
been obtained in this investigation on the syneresis of sodium oleate gels 
in pinene as well as in xylene (cf., Columns 3 and §, Table I). 


Prasad and co-workers*® have found that the relation X”= kt in which 
n and k are constants, holds good in the case of sodium oleate gels in pinene. 
Hence the graphs of log X against log t have been drawn, using the data 
given in Table I, and the plotted points appear to lie on a pair of intersecting 
straight lines (cf., Fig. 2). This indicates that (i) two stages are involved in 
the process of syneresis of gels in xylene and pinene and (ii) the velocity 
of the second stage is much smaller than that of the first. The existence 
of the second stage is noticed only after a sufficient amount of the syneretic 
liquid has been exuded out of the gel. This happens much earlier in xylene 
than in pinene. The existence of the second stage has not been observed 
by Prasad and co-workers® since the period for which they measured the 
syneresis was not sufficiently long to bring out this phenomenon. 
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The above observations can be explained as follows: When a solution 
of sodium oleate in pinene or xylene is cooled, at some temperature the soap 
goes over to the colloidal state due to super-saturation. On further cooling, 
the colloidal particles of the soap agglomerate in the form of fibrils. 
During the formation of colloidal particles and their agglomeration, a part 
of the dispersion medium gets attached to the particles due to solvation. 
This part of the dispersion medium corresponds to the bound liquid in gels 
(cf., Hardy’). The fibrils form the gel structure enclosing the remaining 
dispersion medium (called free liquid by Hardy) within the interfibrillary 
spaces. When the gel starts synerising either due to the collapse of the 
structure as a whole or due to the fibrils coming closer to each other, the 
liquid within the interfibrillary spaces gets exuded; this corresponds to the 
initial stage of rapid syneresis. Subsequently, due to the structural changes 
taking place in the gel, probably a breakdown of the individual fibrils them- 
selves, some of the bound liquid is also exuded; this corresponds to the 
second stage of syneresis. Since the structural changes in a gel take place 
comparatively slowly, the second stage of syneresis is a slow process. 


One of the authors (V. S.) is grateful to the authorities of the University 
of Bombay for awarding him a University Research Studentship during the 
progress of this investigation. 


SUMMARY 


The effect of time on the syneresis of sodium oleate gels in pinene and 
xylene has been investigated. The graphs of synereticum against time are 
smooth rising curves and are not S-shaped as is the case with silica gels. 
The rate of syneresis is large in the initial period and falls off continuously 
with time. The progress of syneresis does not follow the law of first order 
reaction observed in the case of geranin gels. The plots of log X against 
log t for the gels in both the solvents give, in each case, a pair of intersecting 
straight lines indicating the existence of two different processes of syneresis. 
An explanation has been advanced for this behaviour. 
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1. INTRODUCTION 


In Part I, an account of the general theory has been given together with an 
example of the application of the theory to a hypothetical structure. In 
this part, it is proposed to apply the theory to calculate the valucs of rotatory 
power for two crystals of the cubic system, NaClO; and NaBrO;. We shall 
retain the symbols used in Part I. 


Both NaClO,; and NaBrO, belong to the tetartohedral (T) class of the 
cubic system. Their structures have been determined by several workers 
(Dickinson and Goodhue, 1921; Kolkmeijer, Bijeovet and Karsson, 1921; 
Kiby, 1923; Zachariasen, 1929). All of them agree in placing the structure 
in the space-group T*—P2,3 and in giving the number of molecules per unit 
cell as 4. The co-ordinates of the different atoms in one of the enantiomers 
are as follows with respect to a right-system of co-ordinates: 


4Na uuu; 4+4u,4—U,u; usgtut—u; 4—uutt+u 
4C1(4Br) vvv; ete. 
120 xyz; 4+%,4-—y,2; x4+ 94-2; $—x,y,44+2 
yzx; etc. 
zxy; etc. 


The different workers, however, give slightly different values for u, v, x, y, z. 
Of the various measurements, that of Zachariasen is the most accurate; 
but he has studied only NaClO;. Consequently, Zachariasen’s data are 
taken for this crystal, viz., u =0-064, v =0-417, x =0-303, y=0-592, 
z=0-500 and lattice parameter a=6-57A. The data of Dickinson and 
Goodhue for NaClO, agree closely with this and therefore their measure- 
ments for NaBrO, have been used for our calculations. They are: u=0-09, 
y=0-41, x=0-30, y=0-60, z= 0-47 and a=6-71 A. 


The unit cell is shown in Fig. 1, where the four O, groups have been 
labelled 1, 2, 3, 4 respectively. We shall first consider sodium chlorate. 
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2. CALCULATION OF THE POLARISABILITIES 


From the refractive index (np = 1-514) of sodium chlorate, the molar 
refraction of the compound is calculated to be R=12-8. Now, it is well 
known that the positive ion has a low molar refraction (Fajans and Joos, 
1924; Pauling, 1927). For Na, this is of the order of 0-5. As regards 
ClO, ion, no values are available in the literature for the contributions of 
Cl and O atoms separately. However, the following table for ions of the 
type AO, taken from the data in Fajans and Joos’ paper is helpful. The 
table shows that the atom A in the ion AO, always has a small molar 
refraction. 
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Fic. 1. Unit Cell of Sodium Chlorate 
TABLE I 
Ion Ro Ry 
‘CO, | 4°08 0-03 
SiO; | 442 0-10 
NO; - 3-66 0-02 
SO, veel 3°65 0-05 
ClO, «| 3-32 0-04 








We may therefore choose similarly a low value 0-05 for the Cl atom in ClOs. 
Thus, we may take Ry, = 0°45, Rp, = 0-05, Ro = 4:1, whichlead to the value 
12-8 for NaClO,. This set of values is reasonable, since the value 4-1 for 
Ro is about the mean value found in different ions of the type AOs, although 
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larger than in ClO,. Using Heaviside units, we have the polarisability 
a=3R/N, where N is the Avogadro number. Thus, 


ay, = 0°22, Ac, = 0-025, ao = 2:03 x 10-23 (1) 


Now, in calculating the interactions of the different atoms in the structure, 
we may take the ClO; group as a unit. Further, since the Cl atoms contri- 
bute only about 4% of the polarisability of this group, we may neglect them 
for convenience and consider only the O,; group. Calculations from the 
structure show that the O, group is an equilateral triangle, the O—O dis. 
tance (d) being 2-38 A. Consequently, we may calculate the principal 
polarisabilities of this group, using the method of Bragg (1924). If a, 
and «,, are respectively the polarisabilities normal to the plane of the O, 
group and for directions in the plane, then their values are calculated to be 


a, =4-92 x 10°78, a, =7-14 x 10-**. (2) 


In the NaClO, structure, therefore, we may replace each set of 3 oxygen 
atoms by a Single anisotropic group having the above components for the 
polarisability and situated at the centre of gravity of the triangle. The 
co-ordinates of the centres of gravity of the four O, groups, as well as the 
orientation of the normals to their planes, are indicated below: 


Ww, W, W y 4+w,4—w,w 3. w,$+w,4—w 4. +—w,w,4+w 


[111] (11 [li] fi 11] 
with w=0-465. 


We shall make a further simplification by neglecting all the sodium 
atoms also (since these contribute only about 4% to the total polarisability). 
Thus, for the purpose of our calculations, we are left with four anisotropic 
O, groups, | to 4, per unit cell. 


Tt is convenient to take Ox, Oy, Oz along the cubic axes of the crystal. 
Since the O, groups all have their normals along one or the other of the four 
[{111}] axes, their polarisabilities a,; referred to cubic axes can be described 
by two parameters A=(a, + 2a,)/3 and B=(a, —a,)/3. The compo- 
nents of the polarisability tensor for the four O, groups are given below: 


Type J Type 2 Type 3 Type 4 
A B B A-B-B A-B B A B-B 
B A B -B A B —B A-B B A-B_ (4) 
B B A -B B A B-B A -B-B A 


Clearly, A= 6:40 x 10-** and B= — 0-74 x 10-8 (5) 











312 G. N. Ramachandran 


3. CALCULATION OF THE INTERACTION OF NEIGHBOURS 


In the structure, each O, group is surrounded by twelve others (similar 
to the 12 next nearest neighbours of the same type in the NaCl structure). 
The co-ordinates of the twelve neighbours of group 1 are shown in Table I, 
taking that particular group as origin. These are divided into three sets 
of four each, all the groups in one set having the same orientation. Similar 
tables can be prepared for the neighbours of groups 2, 3 and 4 


TABLE II 
Neighbours of O; Group | 





Type Co-ordinates 





2 0:50, —0-43, 0-07; —0-50, —0-43, 0-07; 0-50, 0-57, 0-07; —-0°50, 0-57, 0-07 
3 0-07, 0-50, -0-43; 0-07, —6-50, —0-43; 0-07, 0-50, 0-57; 0-07, —0-50, 0-57 


4 -0-43, 0-07, 0-50; —0-43, 0-07, —0-50; 0-57, 0-07, 0°50; 0-57, 0-07, -0-50 








Now consider light having its electric vector parallel to Ox being 
propagated along the positive z-axis. Then the phase of the wave at an atom 
(x, y, Zz) is exp. (— kaz), since x, y, z are expressed as fractions of the unit 
cell dimension a. Since none of the components of the polarisability tensor 
are zero, all the four types of O; groups will have components of electric 
moment along Ox, Oy and Oz. Taking these into account, the induced 
fields F,, F,, F, are calculated at each group. Each of these three compo- 
nents would produce an additional moment parallel to all three axes. In 
this way, the final induced moments py, and pw, (taking into account the 
influence of neighbours) is calculated for all the four groups in the unit cell. 
Adding them up, we obtain 


Md — . 23__ 2 a 2 mn ‘ 2 
Eye = [25 6 x 10-*— *, (0-115 A?+0-928 AB + 8-50 B )| 

P, 
(E+ 3) © 
Spy = i ms (0-100 AB + 3-096 B?) ke (E+ *) (7) 


‘ 


Since A in the medium is A,/#, where Ay =5-893 x 10° for the sodium D 
line, we have 


P,,=0-9042 (E, + P,/3); P,=0-683 x 10-5” (E, -++ P,,./3) (8) 
From these, we obtain n=1-514 and g/n=9-806 x 10°. Thus 


go 
s 


salle. 
i 8 


= 30°-0 per cm. 
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This compares wel! with the experimentally observed value of 31°-3. The 
excellent agreement between the two must be considered to be fortuitous. 
In fact, small changes in the parameters used would lead to an appreciable 
alteration in the value of the rotatory power. It is however interesting to 
note that the theory gives the rigiit order of magnitude. 


4. ROTATORY POWER OF SODIUM CHLORATE IN THE ULTRA-VIOLET 


Calculations similar to the above can readily be made in the ultra-violet 
if the polarisabilities of the individual atoms are known. Unfortunately, 
all studies on molar refraction, etc., have been confined to a single wave- 
length, viz., the D line and no data are available in the literature for the 
ultra-violet. However, an approximate estimate can be made as shown 
below. We consider the wavelength Ay=2573 A for which the rotatory 
power of NaCiO; has been measured by Rose (1909). The refractive index 
for this wavelength is 1-585, so that the calculated molar refraction is 14-3. 
This has to be divided into the separate contributions of Na, Cl and O. It 
would not be correct to divide them in the same ratio as for A 5893 since 
oxygen would have a smaller dispersion than the other two. A reasonable 
estimate is as follows: 


Ry, =1-0, Rg =0-1, Ro= 4-4. 
Naturally, there is a ccrtain amount of uncertainty in these estimates, but 
they enable us to calculate the order of magnitude of p. These lead to 
ay, = 0:50, ag =0'05 and ag=2-17, x 10-*%, from which we have 
a, = 2-59, a, =1-73, A=6-91, B= — 0-86, x 10-**. 
Substituting in equations (6) and (7) and making the calculations just as for 
\ 5893, one finds p= 226°, as compared with the value 147°-3 found experi- 
mentally. It is interesting to note that the order of magnitude is correct. 
5. CALCULATIONS FOR SODIUM BROMATE 


Rose (Joc. cit.) has also measured the optical rotatory power of NaBrO, 
from 4000 to 7500 A and we give below his data for two wavelengths near 
the upper and lower limits: “i 


A, n p 
7188 1-610 13°-9 
4047 1-643 72°-0 


Taking 7,188 first, the value of Ryap,os is calculated to be 15-0. There appear 
to be no data in the literature as to the individual values of Rg, and Ra in the 
bromate ion. We, therefore, make the following estimates. 


Ry, =0:5, Rg, =0°7, Ro=4°6. 
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Rx, is the same as in NaClO,. Since Br is a much larger atom than Cl, we 
have taken Rg, to be fairly large. Ro has a value larger than in ClO;, but 
we cannot reduce it further without abnormally increasing Rz,. The above 
values lead to a, =2-66 x 10-*% and a, = 1-86 x 10-** for the O, group. 
If we take the co-ordinates of Dickinson and Goodhue (/oc. cit.) for the atoms 
in the unit cell of NaBrO;, the O—O distance is 2-51 A and the centres of 
gravity of the O; groups occur at w, w, w, etc., where w=0-46. This value 
of w is close to the value 0-465 in NaClO,, so that we may take the latter 
here also and no detailed calculations need be made, except to put in the 
appropriate values of the polarisabilities in equations (6) and (7). As in 
NaClO,, here also we neglect the effects of Na and Br, as they contribute 
only about 8% to the total polarisability. It is calculated that A= 7-18 x 10-** 
and B= — 0-80 x 10-** and with these p is found to be equal to 22-1° for 
A 7188. 


Taking 4047, Ry,z-os is 15-65. Obviously, the molar refraction of 
oxygen would vary very little from 47188 to A 4047, since oxygen has a 
dispersion frequency in the remote ultra-violet. Further, most of the dis- 
persion must be due to bromine, since NaBrO, absorbs light completely 
below 4000 A, while NaClO, does so only below 2200A. Hence, we 
estimate the contributions of the different atoms as follows: 


Ry, = 0-6, Rs, = 1 ‘0, Ro= 4:65. 


Thus A=7-21 x 10-**, B= — 0-80 x 10-3 and p is calculated to be 71° 
for A 4047. 


6. DISCUSSION OF THE RESULTS 


The calculated and the experimental values of p for the two crystals 
are Shown in Table III. 








TABLE III 
NaClO,; NabrO, 
iets 
p(th) p(exp.) p (th) p (exp.) 
5893 30°-0 31°-3 7188 22°+1 13°-9 
2573 226° 147°+3 4047 71°-0 72°-0 





The agreement should be considered satisfactory, since the values for 
the polarisabilities are not very accurate. 


The structure of NaClO, can be considered to be derived from the 
NaCl structure by a slight distortion in which both the anions and cations 
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are displaced parallel to the body diagonals. If the O, groups considered 
above occupy exactly the points of a face centred lattice (like the Cl’s in 
NaCl), i.e., w=0-50 in (3), then the rotation vanishes. In fact, it can be 
shown that the ensemble of O; groups then form a structure belonging to 
the space-group T,°—Pa3, which should not exhibit optical activity. It is 
the small deviation from this arrangement that is responsible for the optical 
activity of both sodium chlorate and sodium bromate. 


If the O; groups are taken to be isotropic, then clearly B=0O and the 
right-hand side of equation (7) will be zero, so that p=0. Thus, anisotropy 
of the individual groups is essential for optical activity to be present, as was 
mentioned in Part I. 


We have calculated the anisotropy of the individual groups from theory 
in this paper. This is not quite necessary. In fact, they can be obtained 
from data on the scattering of light by solutions of the crystals. Studies 
of this type have been made for ions like SO, and NO, (Sweitzer, 1927), but 
not for ClO; or BrO, ions. If reliable experimental values are available 
for the anisotropy of the cations, then a more accurate comparison will be 
possible between theory and experiment. 


SUMMARY 


The theory of optical activity outlined in Part I has been applied to the 
cases of NaClO, and NaBrO;. The calculated values for the rotatory power 
agree reasonably well with measurements previously reported. 
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INTRODUCTION 


Doss AND AGARWAL! have studied the discharge of the lead ions at the 
dropping mercury electrode by extending Breyer and Gutman’s? technique 
to different frequencies of the A.C. field ranging from 12-5 to 600 cycles 
per second using the oscillograph for measuring the alternating current. 
Doss and Kalyanasundaram® have studied the discharge of the nickel ions 
by a modified technique. The present paper deals with the application of 
this technique to study the behaviour of lead ions. 


EXPERIMENTAL 


The circuit diagram is given in Fig. 1. The source of A.C. is a B.S.R. 
oscillator capable of giving frequencies ranging from 0-16,000 cycles per 
second. 10 Volts from the oscillator is fed on to a potentiometric arrange- 
ment consisting of resistances R, R,, R2, and R;. This arrangement serves 
to give 45 mv. across R, which is superimposed on the dropping mercury 
electrode along with a D.C. voltage obtained from a Cambridge pH potentio- 
meter (with the galvanometer shorted off) used as a potential divider. The 
A.C. component of the pulsating current is obtained after filtering off the 
D.C. across a condenser. This is measured by means of a two-stage ampli- 
fication system consisting of 1 LD5 valves.* R; is adjusted such that the 
potential drop across this gives the same deflection as that of the potential 
drop across R3. Thus the A.C. current passing through the dropping 
electroce can be calculated from a knowledge of the resistance R;. The 
frequency of the oscillator is maintained constant at the desired value 
making use of the Du Mont 208B oscillograph. . 


The dropping mercury electrodes have the following characteristics :— 


Dropping Dropping 
electrode I electrode II 


(Fig. 3) (Fig. 4) 
m. (mg./sec.) ie i 0-69 0-33 
t in sec. in M KCl-aq; Open circuit ae 3°6 
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The experiment is first performed with ! M KCI containing no other 
dischargeable ions. Hydrogen of high purity is obtained electrolytically 
in the usual manner but with an auxiliary cathode helping in ‘ cleaning up’ 
the residual oxygen diffusion.£ This is passed into the solution for about 
fifteen minutes to remove the traces of oxygen present. The A.C. passing 
through the dropping electrode at the different D.C. potentials are measured. 
The experiments are carried out at different frequencies of A.C. and are 
represented in Fig. 2. 


Figs. 3 and 4 represent the results obtained with 1 M KCI containing 
2 x 10-*M lead ions (in the form of lead nitrate). 


(40) 


A.C. current In ; or 
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D.C. voirs 


Fic. 3. Behaviour of lead ion in pulsating fields 


Fig. 5 represents the polarographic curve of lead obtained with dropping 
electrode I. 
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DISCUSSION 


The A.C. current-D.C. potential curves in the system containing the 
indifferent electrolyte only.— 


An examination of Fig. 2 shows that the A.C. current starts with a 
high value at 0-4 V, falls rapidly upto 0-8 V, remains more or less constant up- 
to 1-6 V and thereafter begins to rise. The first fall of A.C. current in the 
D.C. voltage interval of 0-4 V to 0-8 V coincides approximately with the 
region in which the electrocapillary maximum occurs. It is well known 
that the capacity of the electric double layer on the positive side of the 
electrocapillary maximum is much higher than on the negative side. From 
this point of view, it is of interest to calculate the capacity of the electric 
double layer from the observed alternating current values. There are 
complications arising out of the facts that an average current is measured 
and also that the size of the drop itself is changing with time. The rela- 
tionship between the observed average alternating current and the capacity 
of the electric double layer can be formulated as follows :—The capacitative 
impedance due to the electric double layer at the mercury drop can be put 
as equal to K,ct, where c represents the capacity of the electric double layer 
in farads/sq.cm. and K, ¢ is the area of the mercury drop at time ¢t. If R 
is the external resistance in ohms, the average current as registered in the 
measuring equipment will be 





tmax. 
a 
[Re+ a2 J 
j Pi (ky =a | 
av t 


Making rough calculations, on the basis of the above formula, taking R 
to be negligible we get for the capacities of electric double layer at 0-4V 
and 1-2 V (with reference to N.C.E.) at 50 cycles/sec., the values of 67 and 
26/uF respectively. These values are high as compared with the values of 
49/uF and 22/uF per sq.cm. respectively as given by Philpot.’ These high 
values may be partly due to the fact that the mercury drop is growing all 
the time during the measurement or it may be an artefact of some sort. The 
results however generally indicate the large changes in the capacities of the 
electric double layer near the electrocapillary maximum. The difference 
cannot be due to the fact that the effect of resistive impedance is neglected. 
For, the fact that the current gets nearly doubled as one passes from 25 to 
50 cycles both at 0-8 V and 1:2 V shows that the existing impedance is 
mainly capacitative in character and the resistive part is comparatively 
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negligible. The later rise in the curves beyond 1-6 V is perhaps to be con- 
nected with some reductive process coming into play. 


The A.C. current-D.C. potential curves in the system containing 0-002 M 
lead in MKCI.— 


The observations of Breyer and Gutman that A.C. current maximum 
occurs at nearly the half-wave potential is confirmed by the present system 
as well. It is to be seen that the magnitude of the maximum goes on increas- 
ing with the frequency of A.C. even upto 600 cycles/second. This behaviour 
is different from that of nickel reported earlier® as well as from even that of 
lead itself as measured by Doss and Agarwal? using the oscillograph. 
With nickel, even at a frequency of 100 cycles/second, the maximum has 
disappeared. With lead, in the oscillographic mode of measurement, at 
600 cycles/second, the maximum is of a negligible magnitude. Further, 
the magnitude of the A.C. by this method is much lower than that obtained 
in the present work. 


The above facts can be understood by a detailed consideration of the 
phenomenon. This can be discussed from two points of view, viz., the 
molecular kinetic and secondly, the electrical equivalent of the dropping 
mercury electrode. The molecular kinetic point of view may first be taken 
up. At low D.C. voltages, say, of the order of 0-2 V there is no discharge 
of lead ions and hence the system is behaving like a perfectly polarizable 
electrode. The capacity of the dropping electrode is small, its impedance 
large and hence the resultant A.C. current is small. As the half-wave 
potential is reached, the system becomes imperfectly polarizable due to the 
reduction of the lead ions. There are three important factors to be consi- 
dered in this connection: the electrode process, viz., the discharge of lead 
ions and dissolution of lead from the lead amalgam at the surface of the 
mercury drop, the diffusion processes building up the diffusion layers both 
on the solution side and in the mercury drop and thirdly, the external resist- 
ance of the system. It may be pointed out that the diffusion processes are 
more or less similar in all the systems studied, as the diffusion coefficier. ts are 
oft he same order of magnitude. This would not therefore, bring about material 
differences in the behaviour of the different systems. The main cause of 
the differences in behaviour appears to be the rate of the electrode processes 
which may vary several thousand fold in the different systems. For instat.ce, 
thallium has a rate constant of 1-5 units (moles per sq.cm. per second per 
mol. per c.c.) whereas zinc, with KNO, as indifferent electrolyte, has the 
value of 4:3 x 10-°.6 Apart from the Randles technique which gives an 
elegant way of determining the rate constants, a qualitative idea of the same 
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can be had from the slope of the polarographic curves taken under similar 
conditions. A high slope indicates a fast reaction. When the A.C. is inci- 
dent at the half wave potential, there will be a tendency for the electrode 
process to keep in pace with the pulsating field. The result will be a 
pulsating current. As has already been pointed out,® the A.C. component 
of the current cannot be predicted quantitatively from the polarographic 
curve. This is seen clearly from Fig. 6 wherein the theoretical A.C. current- 


ae 


| | 
+—— +—— 
i | 


| 


| 


AC. CURRENT in JLA 





be) 
0.C- VOLTS 


Fic. 6. Theoretical and experimental curve 


D.C. voltage curve (Curve I, Fig. 6) worked out from the polarographic 
curve (Fig. 5) and the experimental curve at 12-5 cycles (Curve II, Fig. 6) 
are given. The observed A.C. current at the half-wave potential is much 
higher than that predicted on the basis of the polarographic curve. Further, 
this difference is much more at higher frequencies of A.C. The magnitude 
of the fluctuations in the current due to the fluctuating field would depend 
on the rapidity with which the electrode process would respond to the 
fluctuations. With a rapid electrode reaction as in the case of lead, the 
electrode reaction follows easily the fluctuations in the potential and hence 
large maxima are observed even at as high frequencies as 600 cycles/second. 
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Nickel, on the other hand, represents a slow electrode process. The Randles’ 
technique has not yet been applied to a system involving nickel. But the 
rate of electrode process can be judged from the slope of the polarographic 
curve at the half-wave potential. According to the theory of absolute re- 
action rates,’ the rate of the cathodic electrode process should fall exponen- 
tially with the fall in the cathodic potentials. This explains the quick rise 
of the polarographic wave near the half-wave potential. But the actual 
slope should depend on the free energy of activation for the particular 
electrode reaction. The slope would therefore be smaller for slow reactions 
and larger for fast reactions. The low slope of the polarographic wave for 
nickel indicates that it is a slow reaction. With nickel, therefore, the A.C. 
maximum disappears at a much lower frequency of the A.C., viz., 100 
cycles/second, at which the system behaves as if the nickel ions are absent 
and the solution contains only the indifferent electrolyte. 


It is to be noted, however, that the response of the electrode reaction 
is dependent on the resistive impedance of the system as it would affect the 
movement of electrons taking part in the electrode processes. A high 
external resistance would bring down the maximum of the A.C. current 
and would lower the frequency at which the A.C. maximum would dis- 
appear. This phenomenon has perhaps a correspondence with the lowering 
of the slope of the polarographic curve at the half-wave potential with the 
increase in the external resistance. With slow reactions as the external 
resistance would have a greater effect, a pronounced A.C. maximum can 
be obtained only by having very low external resistances. With nickel 
(a slow electrode reaction), for instance, even at a circuit resistance of 600 
ohms, the maximum has disappeared at 100 cycles per second whereas with 
lead with over 1,000 ohms? as circuit resistance the maximum is perceptible 
even at 600 cycles per second. In the present work, the persistance of a 
large maximum with lead even at 600 cycles is due to the very low circuit 
resistance (200 ohms) used in the present mode of measurement. A deétailed 
investigation of the effect of external resistance will be published shortly. 


It is to be pointed out that the difference between the current at the 
A.C. maximum and the corresponding current for the indifferent electrolyte 
at the same D.C. potential (Figs. 2 and 3) is a semi-quantitative measure of 
the amount of electrode process in which the itad ions are involved. It is 
to be seen that this difference is small at low frequencies, increases at higher 
frequencies, reaches a maximum and again becomes small at higher fre- 
quencies. The frequency at which the maximum difference is observed 
appears to be a measure of the rapidity of the electrode reaction if the 
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resistance of the circuit is constant. It is of interest to consider the possi- 
bility of determining the rate constant by these measurements suitably 
modified so as to avoid the complications arising from the changes due to 
the growing drop, the large amplitude of the applied A.C. potential as well 
as the average nature of the A.C. current measured. 


Whatever has been discussed above can be pictured in terms of the 
electrical equivalent. The large rise in the A.C. current at the half-wave 
potential is a result of diminution in the impedance of the dropping electrode 
system. This can be pictured as brought about by two parallel shunts, i.e., a 
resistance and a capacity to the electric double layer at the mercury-solution 
interface. The capacity shunt would have a high impedance for a slow 
reaction as in the case of nickel whereas it would have a low impedance in 
the case of lead. This would completely explain the observed frequency 
dependance. It is of interest to examine how the capacity shunt correspond- 
ing to the electrode reaction would get affected by changes in frequency. 


Grateful thanks are due to Prof. K. S. G. Doss, D.Sc., F.R.Lc., 
F.Inst.P., F.A.Sc., Professor of Chemistry, Indian Institute of Sugar Techno- 
logy, Kanpur, for his valuable guidance and discussions. Thanks are also 
due to Prof. J. M. Saha, m.sc., Director, Indian Institute of Sugar Technology, 
Kanpur, for his kind interest and to the Uttar Pradesh Scientific Research 
Committee for a grant which made this work possible. 


SUMMARY 


The discharge of the lead ions at the dropping mercury electrode has 
been studied by the Breyer and Gutman technique at different frequencies, 
This system has revealed interesting features which are different from the 
systems investigated by others. The results have been interpreted taking 
into account the kinetics of discharge of the lead ion at the electrode. 
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